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SOLUBILITY  0?  NONPOLAR  GASES  IN 
POTASSIUl^i  HYDROXIDE  SOLUTIONS 


By 

Satyendra  Kumar  Shoor 
March,  1968 

Chairman:   R.  D.  Walker,  Jr. 

Major  Department:   Chemical  Engineering 

Measurements  of  the  solubilicy  of  nonpolar  gases  in  aqueous 
potassium  hydroxide  solutions  are  reported  for  a  partial  pressure  of 
one  atmosphere  of  soluce  gas.   The  ten^perature  range  covered  is  25  - 
ICG  C,  and  the  electrolyte  concentration  range  considered  is  0  -  6C  Wt.^ 
KOH.  The  dissolved  gases  studied  are  composed  of  roughly  spherical 
molecules  of  various  sizes  and  include  helium,  argon,  oxygen,  hydrogen, 
sulfur  hexaf luorida,  methane  and  neopentane. 

The  experimental  method  involved  stripping  the  dissolved  gas 
from  solutior.  using  a  carrier  gas  stream,  and  subsequent  analysis  of 
this  stream  by  gas  chromatography.   For  all  the  gases  studied  the 

solubility  falls  off  very  rapidly  with  rising  electrolyte  concen- 

_      .         _o 
tration;  tir.us  ror  oxygen  :Ln  5G  wt.^  KOH  at  2p  C  the  solubility  is 

approximately  1^  of  the  value  for  pure  water.   From  the  experimental 
measurement,  calculated  value-  cf  activity  coefficients,  heats  of  solu- 
tion, energies  of  solution,  and  entropies  of  sjl-tion  are  obtained 

xii 


as  a  runction  of  KOIi  concentration,  temperature  and  solute  gas. 

The  above  experimental  study  has  been  used  to  test  various 
theories  of  gas  solubility  in  electrolytic  solutions.   Electrostatic 
theories  which  make  use  of  th^  concept  of  a  continuous  dielectric 
medium  are  shown  to  be  inadequate  in  explaining  the  salting-ouc  pheno- 
mena.  A  cavity  model  based  on  the  scaled  particle  theory  is  proposed, 
This  model  possesses  several  advantages  over  existing  theories  of  the 
solubility  of  gases  in  electrolytic  solutions  and  is  shown  to  predict 
the  solubility  for  the  systems  studied  in  this  work  fairly  satis- 
factorily up  to  an  electrolyte  concentration  of  50^  by  weight.   The 
salting  phenomena  are  shovm  to  arise  not  because  of  the  ionic  charges 
as  indicated  by  earlier  theories,  but  from  the  presence  of  particles 
(ions)  having  molecular  parameters  (size,  polarizability)  different 
than  those  of  water. 
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CHAPTER  1 
INTRODUCTION 

A  knowledge  of  the  solubilicies  of  gases  in  electrolytic 
solutions  is  of  practical  importance  in  a  variety  of  operations  of 
chemical  engineering  interest.   The  solute  gases  for  which  measure- 
ments are  reported  here  were  chosen  to  represent  a  group  of  nonpolar, 
roughly  spherical  molecules,  having  a  wide  range  of  molecular  diameters, 
and  include  gases  of  direct  interest  for  fuel  cell  operation.   Potassium 
hydroxide  solutions  are  particularly  interesting  because  of  their 
use  in  fuel  calls  and  batteries. 

Previous  studies  have  indicated  that  gases  pass  through  a  thin 
film  of  electrolyte  beiore  reaching  electrocatalytic  sites  of  porous 
fuel  cell  electrodes  (1).   In  hydrogen-oxygen  fuel  cells  the  mass 
transfer  process  at  each  electrode  consists  of  a  diffusional  flux  of 
the  dissolved  reactant  gas  and  a  flux  of  water.   The  reactions  are: 

0^  +  2H^0  +  4e^  40H  (Cathode) 
H,  +  20H-*-  2H,0  +  2e  (Anode) 

In  a  recent  s^^udy  of  electrochemical  kinetics  at  porous  gas-diffusion 
electrodes,  Srinivasan,  'nurwitz  and  Bockris  (2)  have  shown  that  con- 
centration of  the  reactant  gases  near  the  gas-electrolyte  interface  is 
a  controlling  factor  in  obtaining  higher  current  densities.   They 

concluded  that  elimination  of  concentration  polarization  would  permit 

3 

current  dens'-ties  which  are  about  10   times  higher  than  those  obtained 

in  the  presence  of  concen.ration  overpotential.   Under  such  conditions 
a  knowledge  of  the  effect  of  temperature  and  electrolyte  concentr^.tion 


on  the  solubility,  and  diffusiv.cy  of  the  reacting  gases  is  useful  in 
estimating  desirable  operating  conditions. 

In  addition  to  such  pragmatic  considerations,  phase  equili- 
brium studies  over  wide  ranges  of  temperature  and  electrolyte  concen- 
tration are  of  interest  in  zesting  theories  of  gas  .olubilicy  under 
conditions  for  which  the  intermolecular  forces  are  quite  different 
(involving  ions,  dipoles)  than  those  occurring  between  nonpolar 


molecules . 


The  won.  reported  here  formed  part  of  a  larger  program  to  study 
the  thermodynamic  and  transport  properties  of  potassium  and  lithium 
hydroxide  solutions.   Bhatia,  Gubbins  and  Walker  (3)  have  studied 
mutual  diffusion  in  aqueous  potassium  hydroxide  solutions.   Other 
aspects  of  the  program  are  continuing,  and  include  phase  equilibrium 
studies  of  gas- lithium  hydroxide  systems  and  diffusion  of  gases  in 
potassium  and  lithium  hydroxide  solutions. 

For  potassium  hydroxide  solutions  the  salcing-out  effect  is 
particularly  large  for  all  nonpolar  gases,  and  the  gas  solubilities 
to  be  measured  are  very  low,  of  the  order  lo"^  to  lO"^  mole  fraction. 
Because  of  these  low  concentrations,  cor.ventionai  methods  for  solubi- 
lity d>.tcrmination  cannot  be  used.   In  uae  present  investigation,  a  gas 
chromatogrcphic  method  has  been  developed  which  can  be  employed  to 
carry  out  solubility  studies  when  the  amount  of  the  dissolved  gas  is 
very  low.   No  previous  studies  of  gas-electrolyte  systems  seem  to  have 
been  reported  in  which  the  concentrations  are  as  low  as  those  experienced 
ir.  the  present  work.  Zn   addition,  onl,  for  one  or  two  electrolytes 
has  a  systematic  e..:perimental  study  been  made  for  a  series  of  gases  of 


widely  varying  physical  characteristics  over  a  range  of  temperature  and 
electrolyte  concentration.   Such  studies  are  essential  for  developing 
and  testing  theories  of  gas  solubility.   A  further  problem  involved 
in  testing  theories  has  been  the  discrepancies  often  existing  between 
the  results  of  various  investigators. 

In  Chapter  2,  a  brief  discussion  of  the  phase  equilibrium 
relations  which  are  applicable  to  gas-electrolyte  systems  is  included. 
A  review  z£   available  theories  of  solubility  has  also  been  presented 
in  this  chapter   and  this  includes  theories  for  both  gas-liquid  and 
gas-electrolyte  systems.   Previous  attempts  to  account  for  salt  effects 
have  been  along  the  lines  of  the  electrostatic  theories  of  Debye  and 
Huckel.   These  theories  assume  the  solvent  to  be  a  continuous  medium 
of  a  uniform  dielectric  constant,  and  apply  only  to  very  dilute 
electrolytic  solutions. 

A  cavity  model  ba„._d  on  the  scaled  particle  theory  has  been 
previously  shown  to  predict  the  solubility  of  nonpolar  gases  in  both 
nonpolar  and  polar  solvents  with  reasonable  success.   This  theory  has 
been  extended  to  electrolyte  systems  in  the  present  work.   And  the 
model  and  its  development  are  given  at  the  end  of  Chapter  2.   In 
Chapter  5,  results  obtained  by  employing  both  the  electrostatic  theories 
and  scaled  particle  theory  are  co~.pared  and  discussed  in  terms  of 
assumptions  involved  in  each  case. 


CHAPTER  2 
THEORY 

2  . 1 Classical  Thermo dynamics  of  Dissolved  Gases 

Consider  a  nonideal  liquid  solution  in  contact  with  a  vapor 
phase,  at  equilibrium,  for  any  component  i 

u.  =  ,u.  (1) 

where  superscripts  G  and  L  refer  to  gas  and  liquid  phases,  respectively, 
and  [^.    is  the  chemical  potential  (partial  molal  Gibbs  free  energy)  of 
component  i.   We  shall  be  coi»cerned  only  '..'ith  pure  gas  dissolved  in 
aqueous  electrolyte  solutions.   The  gas  phase  will  therefore  contain 
the  solute  gas  and  water  vapor,  while  the  liquid  contains  water,  solute 
gas,  ar.u  the  two  ionic  species  (in  our  case  k  and  OH).   Because  of 
the  electro-neutrality  condition  the  ch>imical  potentials  and  activity 
coefficients  for  the  two  ions  can  be  specified  by  a  single  value  for 
the  dissolved  electroT-te.   We  shall  only  concern  ourselves  with  the 
chemical  potential  of  the  solute  gas.   We  label  water,  solute  gas  and 
electrolyte  as  components  1,  2  and  3  respectively.   For  the  gas  phase 

M?  =  iLi°(T)  +  RTlnf .  (2) 

f . 

—  -*  1  as   ?->  0 

Pi 

The  fugacity  f.  -S  related  to  the  pressure  by  the  following  relation. 


1 


f .  =  o.y.P  O) 

1    1-'  1 


The  fu~acity  coefficient  9.  can  be  obtained  froji  the  virial  equation 
of  state 

ln0.  =  2p2  y.p.  .  -  |  pS  y.y,  C.  .,  -  InZ  (4) 

J     -^      J 

In  the  liquid  phase 

11:   =  /i^(?,T)  -  RTln7.x.  (5) 

The  choice  of  the  standard  state  for  the  liquid  phase  is  arbitrary. 
However  ior  gas- liquid  equilibria  it  is  convenient  to  choose  the 
following  standard  states.   For  water  (1)  we  choose  the  pure  liquid 
at  the  same  temperature  and  pressure  as  the  mixture  (pure  liquid 
standard  ^tate) .   For  the  solute  gas  (2)  choose  pure  solute  in  a 
hypothetical  liquid  state  which  is  defined  by  reference  to  the  behavior 
of  this  component  at  infinite  dilution  and  in  the  absence  of  component 
3;  i.e.,  by  reference  to  Henry's  law  as  defined  below.   We  shall  not  be 
concerned  with  the  activity  coefficient  of  the  electrolyte  (3),  so 
that  specification,  of  its  standard  state  is  unnecessary.   Choosing  the 
above  standard  states  iir.plies: 

Lir.-.7  -^   1   a.-.d    LiaT/^  -»■  1 


2 


X,   — i-  i 


1  -  -  -^2 


■:.  ->  0  (6) 


X3  =0 


Henry's  law  is  defined  as: 


Lin 

x^  -^  0     f^  =  K^x^  (7) 


X3  =  0 


where  K^  is  the  Henry's  law  constant  for  the  gas  dissolved -in  pure 
water.   ?rom  ecuatior.s  (1),  (2)  and  (5),  for  the  dissolved  gas 

The  right  hand  side  of  th-  above  equation  is  independent  of  the  compo- 
sition.  Therefore  using  the  standard  state  as  defined  above 

1^  =  InIC  =  -—  ^ —  (9) 

-^/^■K^  ^   RT         RT  ^^ 

Thus  for  any  concentration 

-2  "  '>'2''2'^2  ^^°^ 

or 

09^2  "  '>'2"'^2^2  ^^^^ 

In  (11)  the  Henry's  constant  depends  on  pressure.   To  explicitly 

accounu  for  this  dependence,  differentiate  equation  (9)  with  respect 

to  pressure  at  constant  temperature  and  composition 

Rewriting  equation  (5)  in  the  form 

,a2(P,T)  =  11^  -   RTln7,X2  (13) 

and  differentiating  with  respect  to  pressure  at  constant  composition 
and  temperature  gives 

*  L 

J  J 


T,n. 


where  V^  is  partial  molal  volume  of  corr.ponent  2. 

Now  consider  a  solution  such  that  :<^  — ►  0  and  x  =  0. 

Equation  (15)  becomes 

where  V   is  partial  molal  volume  of  component  2  in  pure  water  at 
infinite  dilution.   However,  since  the  left  hand  side  of  equation  (16) 
is  independent  of  composition,  it  r.olds  for  any  composition.   Hence 


[   a?    L     ~RT  (^^) 

Integrating,  or.e  cbt^-ins 

?  -0 

^2  =    h'^J    iT'^^.  (18) 

"o 

where  K,  denotes  the  Henry's  law  constant  at  zero  pressure.   Equation 

(11)   can  now  be  ^\rritten  as: 

P-0 

02^2   =  ^2^2^2  ^^P  IJif-  ^0  ^^^^ 

0 

Also  from  equations  (12)  and  (15) 


/Sln^x       V^    /  ^1^2\ 

V  °?    >^  ^   -  RT  -  i    ap    ) 


(20) 


T,n. 
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Using  equation  (10),  equation  (20)  becomes 


T,n, 

J 


This  equation  can  be  used  to  determine  V^  if  experiments  are 
performed  at  different  pressures. 

The  above  equations  are  rigorous.   However,  for  low  pressures 
such  as  those  used  here,  the  treatment  may  be  simplified  somewhat.   At 
1  atm.  pressure  for  nonpolar  gases  the  error  introduced  in  assuming 
0^  =  1  is  negligible.   Also,  at  such  low  pressures  the  pressure 
dependence  of  the  Henry's  law  constant  may  be  ignored,  so  that  the 
integral  in  equation  (19)  may  be  omitted.   Thus  equation  (19)  becomes 

py,  =  72SK2  (22) 

Equation  (22)  may  be  used  to  obtain  7^  for  gas  dissolved  in  an  electro- 
lyte solution  when  K^  and  :-:^  are  known  from  experiment.   IC  may  be 
found  from  experimental  -eta  for  the  solubility  of  the  gas  in  water, 
in  the  limit  of  infinite  dilution  (where  7^  ~  •'•)'  using 

Py2  =  ^2^1  (23) 

Henry's  law  should  apply  at  finite  concentrations  provided  that  the 
concentration  of  dissolved  gas  is  sufficiently  low  that  intermolecular 
forces  between  ^as  molecules  may  be  ignored.   Since  the  molecules  studied 
here  are  nonpolar,  the  intermolecular  forces  between  gas  molecules 
f^ll  off  rapidly  wizh  distance  (as  r   for  dispersion  forces),  so  that 
Henry's  law  should  be  a  good  assumption  for  moderate  gas  partial 


pressures,  where  the  mole  fraction  of  dissolved  gas  is  in  the  range 
lO"^  to  lO"^. 

Partial  Molal  Keat  of  Sol,:tions: 

Differentiating  equation  (9)  wich  respect  to  temperature  at 
constant  pressure  and  composition  yields 

(24) 


V  ^^T    ••p^n_    R  ^  oT   /'  p  ■  R   dT 


Again,  -*  L 

^0    (p  x)  ^'> 
—   '   =  -|  -  Rlir/2  -  Rlnx^  (25) 

and  differentiation  of  equation  (25)  with  respect  to  temperature  at 
constant  pressure  and  composition  leads  to 


/  P,n.  /  P,n 


(26) 


J  J 


H      /  oln7, 


2   .  ^  °""^2  \ 
T      ■         P,n. 

where  H^  is  the  partial  molal  enthalpy  of  component  2  in  the  solution, 
Similar  differentiation  of  equation  (2)  gives 

d.U2(T)/T     .  ^i£/T  V         .  oln0^. 

~^ =  (  ^F-  )     -  (  ->r^J  (28) 

P,n.       ^"  /  P,n. 
J  J 

Assumir.g    the  gas    phase   to   behave    ideally,    we   can  vrrite 


d/^°(T)/T  K^ 


C- 


(29) 
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Since  left  hand  side  of  equation  (29)  is  not  a  function  of  composition, 
H2  remains  unchanged  as  y_  approaches  unity,  and  we  have 

d,u°(T)/T     h? 

'  ■    -If —  =  -  i  ^2°> 

1 

Q 

where  h„  is  the  enthalpy  per  mola  of  the  pure  solute  gas  at  a  pressure 
sufficiently  low  for  the  vapor  phase  act  as  a  perfect  gas, 

^Chen  equations  (24),  (27)  and  (30)  are  combined  one  obtains 


^  SlnK2  ^^        H^    h^     ^  5ln72  \ 

V"^/P.n.  ^■^"?"  ^"^^  l.n. 
J  J 


(31) 


It  may  be  noted  that  both  sides  of  equation  (31)  are  independent  of 
coniposition.   Rearrar.jing  equation  (31),  one  obtains 


/  Sln(7-K^)  V         (::'r  -  h^)     AH 


2    '2'    _  2 


2  2 

P,n.        RT         RT 


(32) 


Using  equations  (10)  and  (32)  yields 


f 

Sin  — 
X2 

> 

(h'^  -  ^2) 

AK^ 

oT 

■ ; 

2 

2 

-    ?,n. 

RT 

RT 

J 

(33) 


where  AH_  is  the  partial  ir.olal  enthalpy  change  on  dissolving  1  mole 
of  ideal  gas  in  the  electrolyte  solution.   If  equation  (27)  is  applied 
to  a  solution  such  that  as  x„  — »•  0,  x  =  0  and  72  ~*  ^  "^  obtain 

r  a/i?/T  s      H 
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where  H^   is  the  partial  molal  enthalpy  of  component  2  in  pure  water 
at  infinite  dilution.   Equation  (32),  then  becomes 


-C         G 

-0 

olnK2 

^2   -   ^2 

AH, 

dT 

'■) 

2 

?,n. 

RT 

RT 

J 

(35) 


where  AH„   is  the  partial  molal  enthalpy  change  on  solution  at  infinite 
dilution  in  pure  x^/ater.   This  equation  holds  generally,  since  che  left- 
hand  side  is  compo sit  ion- independent . 
Partial  Molal  Entropy  of  Solution: 

The  partial  molal  free  energies  of  component  2  in  the  gas  and 
liquid  phases  i-re  related  by 

4  "  ^2  =  (H^  -  H^)  +  TCS^   -  S^  )  =  0  (36) 

Assuming  the  gas  phase  to  be  idei:.l,  we  can  write 

-L     G     (^2  -  <^ 

^2   -  -2  =  -  -H ^  (2^^ 

where  AS^  is  the  change  in  entropy  when  one  mole  of  ideal  gas  dissolves 
in  the  electrolyte  solution. 
Partial  Holal  Energy  of  Solution: 

The  partial  molal  enthalpy  of  component  2  is  given  by: 

^2  "  ^  =  ^^2  -  "^P  ^  P^^^   ■  ^   ^  (3S> 

and  for  an  ideal  gaseous  phase 

^2   "  "4  =  ^^2   "  4^  "  ?^^   "  %    ^  ^29) 
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(e\   -   E^)  -  pV^   (since  V^  «  V^  ) 


AE^  =  AH^  +  RT  (40) 

where  ZiE„  is  the  change  in  internal  energy  when  one  nole  of  ideal  gas 
dissolves  in  the  electrolyte  solution. 
Excess  Functions: 

The  excess  free  energy  of  solution  is  given  by 


where  G"  is  the  free  energy  of  nonideal  mixture; 

/^ro. id  .    ,   .  ^  .  ,   ,   . 

(G  )    IS  the  free  energy  of  ideal  mixture. 


Now     G  =  1  n.u. 
r  1 


=  2n.u'.'(P,T)  +  RT  2  n  liT/.x. 
^  I'^i   '        ^  1   'i  1 

and   (G°)^^  =  Z  n.u'.'(P,T)  +  RT  Z  n.lnx. 
^  I'l  .11 


.   G^  =  RT  2  n . In7 . 

1   1 

The  partial  molal  excess  free  energy  is  given  by: 

3G^ 


G^  =  f  ^  ^        =  RTln7. 
^    ^  °'\  J   T,P,n. 


The   excess    entropy   is 
E 


=     '       jT/  „  =  -  -^   (RTSn.ln-/.)^ 

V  p,n.  oT  i      'i   P,n. 

J  J 

=       .2  n.  irr/.    -   RT2  n.   (   -^^    ) 


P.n. 


The  partial  -lolal   excess   entropy   is 


-E  /  ^^^-i   \ 

;r   =  -Rln7.    -  RT  (       -^   ~     ) 
1  '1  V      oT        ,/. 


f  ,n. 


The  excess  enthalpy  is 

,E   _  E 


H   =  G  -r  TS 
=  -  RT 


^  ^i  i.  -bT-  I 

P.n 


The  partial  nolal  excess  enthalpy  is 


2  Z^^"^- 
Hf  =  -  RT 

1 


2  .2 TV.eories  of  Gas  Solubility 

Theories  of  gas  solubility  in  nonelectro lytic  liquids  have 
been  reviewed  by  Battino  and  Clever  (4).   Recently,  Conway  (5)  has 
sununarized  the  important  modifications  of  the  electrostatic  theories 
for  electrolytic  systerr.s.   The  discussion  of  the  theoretical  rriacerial 
presented  here  is  in  four  sections,  and  reflects  the  different 
approaches . 

The  first  section  deals  v/ich  the  electrostatic  theories.  Most 
of  these  theories  apply  only  to  dilute  electrolytic  solutions,  owing 
to  the  unjustified  assumptions  ir.sde.   The  majority  of  these  theories 
assume  the  solvent  to  be  a  continuous  medium  having  a  uniform  dielectric 
constant.   However,  the  unique  structure  of  water  results  in  a 
characteristic  organization  of  water  r.olecules  (a  hydration  sphere) 
around  an  ion.   In  between  the  hydration  sphere  around  the  ion  and  the 
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main  body  of  the  solvent,  there  may  be  a  transition  zone  which  does 
not  have  any  structure.   These  difficulties  become  more  serious  when 
the  ions  approach  each  other,  i.e.,  in  concentrated  electrolytic 
solutions.   Also  at  higher  concentrations,  the  higher  terms  in  the 
series  expansion  of  a  Boltzmann  exponential  distribution  factor  cannot 
be  neglected.   Nevertheless,  the  theory  is  remarkably  effective  for 
calculating  the  limiting  laws  for  the  deviations  from  ideal  behavior 
as  the  solute  concentration  approaches  zero. 

The  second  section  deals  with  the  regular  solution  theory. 
This  theory  gives  quite  good  results  for  the  solubility  of  gases  in 
nonpolar  liquids.   In  the  case  of  electrolytic  systems,  however,  dis- 
tribution and  orientation  effects  are  not  random.   There  is  a  large 
volume  change  on  mixing  and  the  assumptions  listed  on  page  33  ^^e  '^° 
longer  valid.   In  view  of  the  many  assumptions  needed  in  applying  the 
approach  to  simple  nonpolar  systems,  an  attempt  to  introduce  still 
further  modifications  to  the  theory  so  that  it  may  be  applied  to 
electrolyte  solutions  coes  not  appear  promising. 

The  third  and  fourth  sections  consider  so-called  "cavity  models." 
The  solubility  process  is  assumed  to  be  occurring  in  two  steps.   First, 
a  cavity  of  a  suitable  sizv.  to  accomodate  the  solute  molecule  is 
formed  in  the  liquid,  and  then  the  solute  molecule  is  introduced  into 
this  cavity.   Both  these  steps  involve  energy  changes,   A  cavity  model 
based  on  scaled  particle  theory  is  considered  in  more  detail  in 
Section  2.2.4  and  this  model  has  been  extended  to  aqueous  electrolytic 
solutions.   This  section  contains  the  theoretical  contribution  of  the 
present  work.   The  theory  developed  there  is  tested  in  Chapter  5. 
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2.2,1   Electrostatic  Theories 

The  addition  of  an  electrolyte  to  an  aqueous  solution  contain- 
ing a  neutral  molecule  frequently  tends  to  decrease  the  solubility  of 
the  neutral  species.   This  corresponds  to  an  increase  in  activity 
coefficient  of  the  molecule.   "Salting -out"  of  this  kind  is  not  general, 
and  examples  of  the  reverse  behavior  of  "salting -in"  may  be  found. 
The  electrostatic  theories  discussed  here  attempt  to  explain  these 
phenomena  on  the  basis  of  coulombic  forces. 

A  comprehensive  review  of  electrostatic  theories  in  relation 
to  salt  effects  was  made  by  Long  and  McDevit  (6).   A  review  by  Conway 
(5)  summarizes  the  recent  important  improvements  made  in  these  theories. 
Debye  and  McAulay  (7) 

The  earliest  electrostatic  theory  of  salt  effects  v;as  that 
of  Debye  and  McAulay  (7).   These  authors  made  the  following  assumptions 
as  to  the  nature  of  the  electrolytic  solution: 

1.  The  solvent  miolecules  are  treated  as  a  continuous  die- 
lectric medium  whose  dielectric  constant  is  uniform  throughout  rhe 
solution. 

2.  There  is  a  continuous  charge  density  function  around  a 
given  ion  due  to  the  ionic  atmosphere. 

3.  The  solution  is  very  dilute  with  respect  to  electrolyte. 
These  authors  compare  two  solutions,  each  containing  water, 

electrolyte  and  dissolved  gas,  and  each  having  the  sam.e  total  number 
of  moles  of  electrolyte  and  dissolved  gas.   However,  the  first  solution 
is  infinitely  dilute  (i.e.,  contains  an  infinite  number  of  moles  of 
water),  whereas  the  second  solution  is  the  system  of  interest,  of 
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finite  concentration,  containing  n  ,  n^ ,  ....  n n  ions  per 

cubic  centimeter  of  charges  e, ,  e^  ....  e e  ,  etc. , resulting 

froi?.  n'  molecules  of  electrolyte.   The  first  solution  is  converted  to 
one  having  the  same  intensive  properties  as  the  second  solution  of 
finite  concentration  by  the  following  two  steps. 

1.  Removing  the  charges  e.  on  the  ions  infinitely  slowly. 

2.  Compressing  the  solution  isotherm.ally  and  infinitely 
slowly  to  the  final  volume,  by  means  of  a  piston  semi- permeable  to 
water  and  finally  restoring  the  charges  infinitely  slowly.   The 
electrical  work  involved  in  discharging  the  ions  at  high  dilution  and 
recharging  them  in  the  solution  of  finite  concentration  is  then  con- 
sidered.  Since  the  potential  at  the  surface  of  a  spherical  ion  of 

e . 

radius  b.  in  a  medium  of  dielectric  constant  D  is  -r~—    ,  the  work  of 
J  Db. 

J 
charging  by  a  Debye  charging  process  (8,9)  is 

4=1   2,.,„        2 
s    ^  n    e.c,d^  s  n.e. 

W  =  Z  n.  /  -^  =  Z       JJ.  (41) 

c    .    J  /   Db .     ,   2Do . 

J      e=o       J        ^    J 

where  4  is  a  charging  param.eter  which  varies  from  0  to  1  as  the  charge 
e  increases  from  0  to  e.   Similarly  the  work  of  discharging  is  obtained 

as  2 

s   n.e. 

^d  =  -  ^  2Ff  ^''^ 

J    o  J 

Next,  the  work  due  to  the  ionic  atmosphere  is  calculated.   The  poten- 
tial ip.    at  the  position  of  the  ion  j  due  to  the  ionic  atmosphere, 
obtained  from  Debye-Huckel  theory  (0)  is 


e./c 


-     ^  (43) 
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where  \j2  \ll 

2  \  '   /   11. 


.^  DkT  y    l^  V-  J 


.-C  =    ^T^       2  ^  Z  1  (44) 


where   e.  =  Z.e 

J    J 


Z.  =  valence  of  ion  j 

J 

£  =  electronic  charge 

V  =  volume  of  the  solution 


n.  =  number  of  ions  of  type  1  in  volume  V, 

J 


Work  due  to  the  ionic  atmosphere  is 

W.  =  2n   /— ^ de=§^  2ne  .:  (45) 

1  'e=o  ^ 

The  appearance  of  /c^  in  the  above  expression  results  from  the  fact  that 
K   is  proportional  to  the  charge  e.   As  a  first  approximation  the  total 
electrical  work  obtained  by  the  above  processes  is  set  equal  to  the 

increase  in  free  energy  of  the  system.   Thus 

2  2 

,sn.e.     ,    sn.e.     ,^   s 

i     J        o  1      J  1   J  J 

The  dielectric  constant,  D,  of  the  solution  is  represented  by  the 
following  relation 

D  =  D  (1  -  Bn  -  S'n')  (47) 

o     '     ' 

where  n  and  n'  are  the  numbers  of  molecules  of  nonelectrolyte  and 

electrolyte  respectively  and  j3  and  /3'are  empirical  constants.   D  is 

dielectric  constant  of  the  solvent.   By  substituting  equation  (47) 

equation  (46)  and  neglecting  the  higher  powers  of  n  and  n',  the 
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following   equation   for    free   energy   change    is   obtained: 

2  2^ 

■7        sn.e.  — ,     ,sn.e.  s  2 

^^        2D     f     b.  2D  f  ^r~        3D       f     J    J  ("^S^ 

olj  o-'-J  ol 

Computing  the  contribution  to  the  chemical  potential  of  the  non- 
electrolyte  by  the  addition  of  electrolyte  and  using  its  thermodyna- 
mic relationship  to  solubility,  gives  the  final  equation  of  Debye  and 
McAulay  (7)  which  is 

S°     -    s   n  e  ^ 
z        o   1      J 

v;here  S^  and  S^  are  solubilities  of  nonelectrolyte  in  the  absence 
and  presence  of  electrolyte. 

Randall  and  Failey  (10)  tested  the  validity  of  equation  (49) 
for  several  nonelectrolyte-wauer-electrolyte  systems  up  to  an 
electrolyte  concentration  of  4  molal.   It  was  found  that  concentration 
dependence  of  salt  effects  is  predicted  within  the  limits  of  probable 
experimental  error  in  many  cases. 
Debve  (11) 

Debye  (11)  later  developed  a  more  exact  theory  to  obtain  an 
expression  for  the  concentration  of  a  nonelectrolyte  molecule  in  the 
presence  of  a  spherical  ion  of  radius  b .  as  a  function  of  its  distance 
from  the  ion.   He  calculated  the  free  energy  for  the  real  solution  by 
taking  into  account  the  energy  of  the  electric  field  surrounding  the  ion. 

If  X  is  the  field  strength  in  a  given  volume  element,  the 
energy  per  unit  volume  is  given  from  electrostatics  by 

.2 


(DXV 
SrrD 


(50) 
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where  D  =  dielectric  constant  of  the  medium. 
For  an  ion 

X  =  -  -^  (51) 

Dr 

Thus  the  energy  of  the  volume  element  at  a  distance  r  from  the  ion  is 

given  by 

2 

^    dV  (52) 


8nr^D 


The  total  free  energy  of  the  real  solution  is  obtained  as 

2 


=  f 


soln 


n^(0^  +  kTlnx^)  +  n^<i0^   +  kTlnx^)  +  ^  ^   i  dV 

Stx   D  -J 


(53) 


where  subscripts  1  and  2  denote  solvent  and  nonelectrolyte,  respectively 
n.  =  number  of  molecules  per  cm.   of  component  i 
0^  =  mJ^/N  (54) 

If  6n^  and  6n„  are  small  variations  in  n  and  n„ ,  then 

6(|>  =   /   I  &nA0.    +   kTlnx^  -  ^-y  -r  ^  )  +  6n_(0-  +  kTlnx^ 
I       i   i        i   „„  '4-  _-i  on,       12.  2. 

^   L  STTr   D    1 


'   \§?-)ldV  (55) 


87^"^  D^  ^^2  ■  J 

By  using  the  condition  of  equilibrium,  645  =  0  or  $  is  a 
minimum  subject  to  the  condition  chat  the  total  number  of  molecules 
does  not  change,  and  assuming  no  volume  change  on  mixing,  the  following 
eouation  is  obtained: 
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v„ln  — r-  -  V  in  — r  =  v  t- -  v  ^ —   —  (56) 

2   ^0     1    0   o^i-n^-n  V   2  on     1  dn   /  4 

where  v.  is  volume  per  molecule  of  component  i,  and  x  and  x  are  the 
values  of  x  and  x„  at  a  large  distance  r.   Debye  (11)  then  treated 
the  salting-out  of  component  2  by  considering 

-|  <  1  and  ^  (57) 

very  large  so  that  equation  (56)  becomes 

0  -(R/r)^ 
X2  =  X2  e 

where  ^ 


(58) 


If  we  assume  a  solution  saturated  with  component  2  (at  low  concentration) 
in  the  presence  of  2  n.  ions  per  unit  volume,  the  average  number  of 
nonelectrolyte  molecules  surrounding  an  ion,  can  be  obtained  from 

X       x°         4 

z     =  r   dV    =    /"^  dV  =  ^  [ ^<^l^">   dV    (60) 
2   J    2  /  v^      v^j 

V  V  V 

The  integration  is  over  the  volume  V,    at  the  disposal  of  an  ion. 
Similarly,  the  number  of  molecules  of  component  2  in  the  volume 
element   before  the  addition  of  electrolyte  (which  is  the  same  as 
the  value  of  Z„  at  very  large  r)  is: 


Z 


0 
P  _  ^2 
2  ~  V,  / 


fdv  (61) 


Using  equations  ( 60)  and  (6l)  and  recalling  that  2  n.  /  dV  =  1,   (62) 
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the  following  equation,  for  salting -out  is  obtained: 

^=1   -   Zn.      Ul   -   e-(^/^^)dV  (63) 

^2  V 

For  dilute  solutions,  this  becomes 

s  ""  4 

-|=  1  -  2n.  ^'  (1  -  e"^^/^^  )1.7:r2dr  (60 

^2  b. 

J 

Debye  (ll)  simplified  equation  (64)  by  expressing  dielectric  constant  D 
as  a  linear  function  of  the  mola  fractions.   The  integral  involved  in 
the  equation  (64)  was  evaluated  by  Gross  (12)  to  give  a  series  solution. 

Long  and  McDevit  (6)  compared  experimental  activity  coefficients 
with  those  calculated  from  Debye 's  theory  for  a  variety  of  nonelectro- 
lytes  in  aqueous  electrolyte  solutions  at  electrolyte  concentrations  up 
to  1  molal.   Their  conclusion  vjas  that  the  theory  gives  results  of  the 
correct  order  of  magnitude.   However,  the  theory  cannot  reasonably  account 
for  the  marked  variation  in  the  effect  of  different  electrolytes.   Long 
and  McDevit  (6)  attributed  discrepancies  between  the  experimental  and 
theoretical  values  to  the  fact  that  Debye 's  theory  considers  only  the 
primary  electrical  effects  and  does  not  include  the  "displacement"  and 
"structural"  contributions.   Long  and  McDevit 's  argument  is  supported  by 
the  fact  that  the  theory  usually  gives  the  closest  predictions  for  just 
the  cases  where  the  latter  factors  should  be  minimized,  e.g.,  sodium 
and  potassium  chlorides. 

Gross  (13)  extended  the  theory  of  Debye  (ll)  to  include  the 
effect  of  the  ionic  atmosphere.   This  effect  requires  that  the  energy 
of  the  volume  element  at  a  distance  r  from  the  ion  increases  by  an 
amount 


2  r  ■\      /     -  '^r 
e 


G     I    e ^j 


SttD  i_   ^r  V,     r 
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2  2,,     ,      ,  ,    -iKr 

dV   =   ^   .^^    Y"^^ dV  (65) 

Sttt  D 


where  k;  is  a  parameter  defined  by  equation  (44)  and  b  is  the  ionic 

radius . 
Butler  {\K) 

Butler  (14)  obtained  a  relation  very  similar  to  that  of  Debye 

and  McAulay  (7)  by  a  slightly  different  procedure.   However,  it  is 

easier  to  use  since  it  is  expressed  in  terms  of  more  accessible 

molecular  quantities.   He  calculated  the  electrical  work  required  to 

bring  a  molecule  of  component  2  from  infinity  to  a  distance  r  from  the 

ion  by  the  expression 

W  =   /  P„6v.dX  (66) 

J   2  2 

where  P^  =  polarization  produced  per  unit  =  a  X 

a^  =  polarizability  of  molecule  2 

X  =  field  strength 

6v^  =  volume  of  molecule  2 

Thus,  the  total  electrical  work  to  bring  a  molecule  of  component  2 
from  infinity  to  a  distance  r  from  the  ion  and  displacing  an  equal 
volume  of  component  1  away  from  the  ion  is 

(a,  -  a„)e 
W  =  ^, 5v-  (67) 

Next,  the  Boltzmann  equation  for  the  distribution  of  molecules  of 
nonelectrolyte  was  assumed  to  obtain 


^2  =  \   ^^H  — E 7JT  ^^2 

L  2D  r     . 


(68) 
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Frora  equation  (68)    the  amount  of  nonelectrolyte  salted  out  by  a 

3  0 

single  ion  per  cm.,  T     is  obtained  by  integrating  (Z^  -  Z^),  the 

3 
deficit  in  the  number  of  molecules  of  component  2  per  cm  .  at  a  distance 

r,  over  the  total  volume.   Thus 

<r:  2D  r 

b 

In  the  above,  only  linear  terms  in  the  expansion  of  the  exponen- 
tial of  equation  (68)  have  been  used — an  approximation  which  is  justi- 
fied only  at  very  low  electrolyte  concentrations.   For  the  solution 
containing  m  g . ions  of  charge  e  ,  m^  g-ions  of  charge  e^    ....per  liter 
equation  (69)  becomes 

0  2  ~ 

S„  -  S„   a,  -  a„  27r6v    e-  m-N 

_2 £  =  _i £ ^  Z  — ^^ (70) 

0        kT      2     b.  1000  ^'^ 

S^  D    J   J 

Belton  (15) 

Belton  (15)  modified  the  equation  of  Debye  and  McAulay  (7)  for 
the  case  of  a  polar  solute  gas  in  a  polar  solvent.   He  added  the  work 
done  due  to  the  permanent  dipole  moments  of  the  polar  components  to 
the  electrical  work  obtained  in  the  case  of  nonpolar  components  (14), 
Kis  final  equation  for  solubility  is 

S°   471  a^o^n.  47Tb2  Z.n.e 

^°S  F  =  —7^ ^"-l   -   ^2>  ^  ~XkT^  ^^1  ■  ^2^ 

2    2D^bkT  0 

where  b  is  the  distance  between  the  centers  of  the  ion  and       (71) 

the  nonelectrolyte  molecule  at  the  position  of  the  nearest  approach,  and 

'fi^    and  ,a„  are  the  dipole  moments  of  solvent  and  nonelectrolyte  solute. 
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respectively.   The  author  used  equation  (71)  to  calculate  the  activity 
coefficient  of  five  solutes  in  different  aqueous  salt  solutions,  but 
did  not  obtain  good  agreement  with  the  experimental  values. 
Bockris  et  al.  (16) 

Bockris  _et  al .  (16)  sought  an  explanation  of  the  salting-in 
effects  for  benzoic  acid  by  tetra  alkyl  ammonium  iodides  in  aqueous 
and  ethylene  glycol  solutions  by  considering  the  dispersion  forces 
between  the  ion  and  the  nonelectrolyte  molecule  in  addition  to  the 
coulombic  Jorces  considered  by  previous  workers.   They  obtained  the 
following  expression  for  the  work  done  in  bringing  1  molecule  of 
nonelectrolyte  from  infinity  to  a  distance  r  from  the  ion  and  dis- 
placing an  equal  volume  of  solvent  away  from  the  ion: 

v^        9  9        ^2 

(a  ^  -  a  )z  e        (^^  —  "  ^{) 

1  1 

-"■—7^ * 1 ('^> 

2D  r  r 

where  subscripts  1  and  2  denote  solvent  and  nonelectrolyte  solute, 
respectively,  and  \.    is  a  constant  for  dispersion  forces  between 
component  i  and  the  ion,  obtained  from  London's  dispersion  formula. 

Using  equation  (72)  and  following  the  procedure  of  Butler  (14), 
the  following  equation  for  the  salt  effect  is  obtained: 

^2         2~ 
„    „0        (a, a^)27TG  N 

^2  -  ^2         r  _ill__l— _  1 

$2         i-    1000  D^kT      J 


m. 

-^  (73) 

b-r 

1 


where  m.  is  g.  ion  of  charge  Z.£  per  liter  of  solution, 
1  X 


„2 

^  .m. 
1    1 

r           V 

Attn        -| 
-  X^)    3000kT 

b. 

X 
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These  authors  compared  experimental  values  of  (S„  -  S^)/S„rn 
with  those  calculated  from  equation  (73) for  salting-in  of  benzoic 
acid  by  five  different  salts  in  aqueous  solutions;  they  found  that  the 
theory  gave  the  correct  sign  and  order  of  magnitude  for  the  various 
systems.   However,  they  pointed  out  that  equation  (73)  will  not  hold 
good  at  high  concentrations  because  of  the  unjustified  expansion  of 
the  exponential  distribution  factor  and  the  neglect  of  the  ionic 
atmosphere. 
Altshuller  and  Everson  (17) 

Altshuller  and  Everson  (17)  modified  the  theory  of  Debye  and 

McAulay  (7)  by  taking  into  account  the  variation  of  the  parameter  /c 

v/ith  concentration  of  nonelectrolyte.   Their  final  equation  is 

0   —   2~  2 

S     ,£D  £  N     ^  s  m.Z.       s        N 

In  ^  =  -^ 5—  (  2  -i-^  -  K  Z    m  z/  )  (74) 

2   2000  D~kT    1     i      1 

These  authors  compared  the  experimental  results  for  x^ater-ethyl  acetate- 
alkali  halide  systems  up  to  electrolyte  concentrations  of  1  molar 
with  those  calculated  from  equation  (74),  and  with  the  equations  of 
Debye  and  McAulay  (7)  and  Debye  (11).   It  was  found  that  none  of  the 
equations  were  fully  satisfactory.   Altshuller  and  Everson  suggested 
that  hydration  radii  and  not  the  crystallographic  radii  should  be  used 
in  the  equation  of  Debye  and  McAulay.   This  would  remove  the  objection 
to  the  latter  equation  that  electrical  work  tends  to  infinity  where 
crystallographic  radii  approach  zero. 
Mikhailov  (18) 

Mikhailov  (13)  introduced  the  concept  of  a  salting  out  sphere. 
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the  radius  ot  which  is  defined  as  follows: 


where  Si.  is  the  number  of  ions  per  unit  volume  of  solution.   He 
1 

pointed  out  that  an  allowance  should  be  made  in  the  derivation  of 
salting-out  equations  for  the  reduction  in  the  salting-out  sphere  as 
the  electrolyte  concentration  increases.   This  implies  that  the  integra- 
tion in  equations  (64)  and  (59)  should  be  carried  out  from  b.  to  r  , 
rather  than  from  b  to  oo,  which  is  justified  only  at  infinite  dilution. 
Thus, 

Z.  =   /  Z  e  "■    477r-dr  (76) 

'b. 

1 

where  W.  =  work  term  given  by  equation  (72). 

Mikhailov  \>rrote  the  work  in  the  form 

E    D 
W.  =-i+-|  (77) 

r    r 

where  E.  and  D.  are  given  by  the  first  and  second  term,  respectively, 
of  equation  (72)  without  r.   The  integral  in  equation  (76)  is  evaluated 
by  dividing  it  into  two  parts;  from  b.  to  oo,  and  from  r   to  oo,   a  more 
justifiable  expansion  of  the  exponential  can  then  be  made  for  r  >  r 
(W  decreases  rapidly  for  large  values  of  r) .   In  this  way  the  equation 
for  salting-out  by  1: 1  electrolytes  is  obtained 


2? 


S„       47TNC     ,      ,       p        -W./kT 


2 


2  b. 

1 

~  s  A/3  ^   ~  V  2 

8-N  \-'^  „  4/3   .  /  4-^  V  ,^   ,  ^  ,   2 


1000  ;      ^s       -  M,  3o3o  j    ^\  ^  ^2^  ^3 

(78) 

where  S^  =  E^  =  E 

and    C  =  molar  concentration  of  salt. 
s 

Mikhailov  (18)  used  equation  (78)  to  calculate  the  salting-out  of 
tributyl  sulphate  and  benzene  by  various  salts  up  to  1.5  molar  salt 
concentrations.   The  agreement  with  experimental  results  was  not  good. 
However,  correct  estimates  of  sign  and  magnitude  of  the  deviations  from 
the  limiting  law  were  obtained. 
Conway,  ez   al.  (19) 

Conway,  et  al.  (19)  showed  that  limitations  in  previous 
electrostatic  theories  arise  from  an  unjustified  expansion  of  the 
exponential  distribution  factor  and  the  neglect  of  solvent  dielectric 
saturation  effects  in  the  vicinity  of  the  ion,   where,  relatively,  the 
maximum  salting-out  is  developed.   They  followed  the  method  of  Debye 
(11)  and  made  important  modifications  to  correct  for  the  above  limita- 
tions.  They  started  with  a  distribution  equation  very  similar  to 

equation  (60)  which  can  be  restated  as 

0 

^9  2 

dZ^  =  -=■   exp[-Au/kT]4Trr  dr  (79) 

where  dZ^  =  number  of  nonelectrolyte  molecules  in  a  spherical  volume 
element  of  thickness  dr  at  a  distance  r  from  the  ion. 

Au  can  be  evaluated  by  Debye ' s  (11)  treatment  and  is  given  by 
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To  obtain  the  nuir.ber  of  nonelectrolyte  molecules  in  a  volume  associated 
v/ith  one  ion,  the  expression  must  be  integrated  from  b  to  R,  where  b 
is  the  radius  of  the  ion  and  R  is  a  critical  radius  (very  similar  to 
the  salting-out  sphere  radius  of  Mikhailov  (18))   corresponding  to  the 
volume  available  per  ion  in  the  solution  and  defined  by 

1/3 
R  =  (fr^  '  (81) 


y47lNm  J 


where  m  is  g.  ion  of  electrolyte  per  liter  of  solution.   Thus, 

Using  the  above  equations,  and  following  the  procedure  of  Debye  (11), 
one  obtains  an  equation  for  the  salt  effect  in  the  following  form: 

7°  -  Z     S°  -  S        ~    ^ 

12 tl       11 _2    At^'i^  r,.  -Au/kT,  2,  .... 

2       ^2  b 

These  authors  pointed  out  that  the  exponential  Au  should  not  be 
expanded  to  its  linear  form,  as  done  by  Butler  (14)  and  Bockris  (16), 
because  Au  is  usually  not  small  in  comparison  with  kT  near  the  ions, 
where  the  largest  contribution  to  the  salting-out  effect  is  expected. 
It  vjas  shown  that  there  are  two  regions  in  the  vicinity  of  an  ion.   In 
the  primary  hydration  shell  only  field-oriented  water  molecules  are 
preseni:,  so  that  the  dielectric  constant  of  the  solution  is  very  small. 
Consequently  Au/kT  is  large  and  expL-Au/kT]  is  small  compared  to  unity. 
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Beyond  the  primary  hydration  shell  Au/kT  is  small  compared  to  unity 
owing  to  the  rapid  rise  of  D  with  increasing  distance  from  the  ion. 
Hence  the  exponential  in  equation  (83)  may  be  expanded  retaining  only 
the  first  tern  of  the  series.   Thus,  the  integration  of  the  salting- 
out  equation  is  conveniently  carried  out  over  two  integrands;  from  b 

to  -,  ,  where  D  «  D^,  and  from  r  to  R,  where  D^  D    The  equation  (83) 
n  U  n 

then  becomes 

qO   <;     -  -  ''^        }  n 

!2_lI=itgL.  Tr^dr^  r(Au/kT)r2dr  '  (84) 

^2^  b        r, 

^  n 

To  evaluate  Au,  a  limiting  case  of  Kirkwood's  theory  of  dielectrics 
(20),  applicable  only  at  low  field  strengths,  is  used.   This  is  per- 
missible because  the  field  strength  is  not  very  high  in  the  region 
r  <  r  <  R.  According  to  this  theory 

h 

iooo(D  -  ::) 

^ =  V^Dq  -  (9/2?2  (8d) 

and 

P^  =  4/3  rka^   +  ^  )  (S6) 

where  V  is  the  molar  or  partial  molar  volum.e  of  nonelectrolyte,  and 

ji  is  the  dipole  moment  of  the  individual  nonelectrolyte  m.olecules  in 
solution.   The  parameter  Jl   is  the  vector  sum  of  the  moment  of  the 
central  molecule  plus  that  of  all  the  neighbors  which  may  affect  its 
orientation.   Thus,  for  a  non-associated  polar  molecule,  we  have 

^  =  ^  "^  ,  (V  D,  -  (9/2)?,)  \  ^h^  '^  <  ^  ^^^^ 

'^^   BTikTADQ-     ""  "  r 
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and  equation  (84)  becomes 


'2-'   ^-^      .      3       .3.    zV 


.0      3000  ^^h  "  ^^^  + T  (V,Dp  -  ^/2)PJ(  -^  -  I) 

^2™  2000kTDQ    2  0        2'^^   R^ 

The  first  term  in  equation  (88)  is  proportional  to  the  nulb!r 
of  water  molecules  held  in  the  primary  hydration  shell.   This  term 
is  obtained  by  a  different  but  equivalent  method,  i.e.,  by  considering 
the  salting-out  by  ions  as  being  essentially  determined  by  the  extent 
of  primary  hydration  of  the  ion  insofar  as  the  water  in  the  primary 
hydration  shell  is  effectively  removed  from  its  normal  solvent  role. 
The  first  term  of  equation  (88)  can  then  be  replaced  by 

•  ISnm 


lOOOd  -  mW 
where  n  is  the  hydration  number 

m  is  molarity 

d  is  density  of  solution 

W  is  molecular  weight  of  dissolved  salt 
With  this  substitution  equation  (88)  becomes 


,0       -  lOood  -  „„  -  ^^^^;^2  ("2%  -('/2)P2)(  r-  -  I  ) 

^  0 


(89) 
Conway  et  a^^  (19)  used  equation  (89)  to  calculate  the  salt 
effects  on  rare  gases  and  other  neutral  molecules  by  the  alkali  halides. 
and  they  compared  the  results  with  experimental  values.   It  was  found 
that  the  dependence  of  salt  effects  on  the  size  of  the  nonelectrolyte 
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was  predicted  fairly  well,  whereas  the  dependence  on  the  nature  of  the 
ion  for  a  given  nonelectrolyte  was  predicted  well  only  when  the  non- 
electrolyte  molecule  was  small.   For  large  nonelectrolyte  molecules  the 
agreement  was  poor. 

Recently,  Ruetschi  and  Amile  (21)  introduced  a  procedure  which 
was  very  similar  to  that  used  by  Conway  _et_  al .  (19)  to  obtain  the 
following  equation  for  salt  effects: 


5°  -  S^ 

=  2n.  f  -r  Tn:"  )  +  2  n.   /  (  1  -  exp 


1 y  3    hy         1 


rN'Tn,    2 


KfnT  dr 


(90) 


where  Si.  is  number  of  ion  per  cm. 
1 

r  has  the  same  meaning  as  R  in  equation  (81). 
The  first  term  of  above  equation  is  the  volume  of  the  hydrated  ion 
and  is  evaluated  from  hydration  theory.   .According  to  this  theory  the 
ratio  of  the  solubilities  in  the  pure  solvent  and  in  the  solution  of 
an  electrolyte  may  be  set  equal  to  the  ratio  of  the  volumes  of  the  "free" 
water  without  regard  to  electrostatic  forces.   Thus, 


and 


4            r 

Lim 

n.  -*■  0 

1 

S„        1-2  n.V. 
2                        11 

In  -T-  =  2  n.V. 
^2 

V.  is  the  volume  of  the  hydrated  ion. 
2.2.2   Reg:ular  Solution  Theory: 

Regular  solution  theory  has  been  frequently  applied  to  gas 
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solubilities.   Eildebrar.d  and  Scott  (22)  define  a  regular  solution  as 
one  in  which  orienting  and  chemical  effects  are  absent,  and  in  which 
distribution  and  orientation  of  molecules  are  random  as  in  an  ideal 
solution.   Thus  a  regular  solution  is  one  involving  no  entropy  change 
when  a  small  amount  of  one  of  its  components  is  transferred  to  it  from 
an  ideal  solution  of  the  same  composition,  the  total  volume  remaining 

unchanged.   With  this  concept  of  a  regular  solution,  the  free  energy  of 

M 
mixing,  AF'  can  be  written  as 

M     M      M 
AG   =  AH  -  TAS 

M      M      M 
=  AE  +  PAV  -  TAS 

=  AE^'  -  TAS^  (91) 

where  the  term  PAv   is  assumed  to  be  negligible. 

Hildebrand  and  Scott  (22)  derived  the  following  expression  for  the 

free  energy  of  mixing  of  N^  moles  of  component  1  and  N^  moles  of 

component  2: 

,2 


AE-  =  2TrN^(N^V^  +N2V2)  ,  -r     I   ^n^°ll  -  \ 

1 


+ 


^2  1 


To^^^r^dr  -Vxf^2Si^'dr) 


vl      -■   ^"-       '1''2 

where  subscripts  1  and  2  denote  component  1  and  2,  respectively. 

N.V. 
1  1 
ijj     =■   volume  fraction  or  component  1:   2~n~V~ 
■'■  .   i  i 

K  =  Number  of  moles  of  component  i 

1 


(92) 
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(J).  .  =  potential  betxi^een  the  central  molecule  i  and  each  surround- 
mg  molecule  j 

g..  =  radial  distribution  function  (probability  of  finding 
molecule  j  around  molecule  i  in  the  mixture) 

r  =  variable  distance  from  the  central  molecule 

g.  =  radial  distribution  function  for  pure  component  i 


g-i  =  gii/'i/'i 


The  above  expression  has  been  further  simplified  by  making 
the  following  assumptions: 

1.  g^i  =  gi 

2.  Additivity  of  the  energies  of  molecular  pairs 

■^'    ii    n    6  (93) 

r    r 

where  j  and  k  are  constants  for  repulsive  and  attractive 
forces  respectively 

4.   All  radial  distribution  functions  can  be  expressed  as  the 

]^ 
same  universal  function  g(y)  where  y  =  -r-   and  d   is  the 

d       m 

m 

position  of  the  first  maximum  of  the  radial  distribution 
function  g(r)  in  the  plot  g(r)  vs.  r  (Figure  l) 

^-    '^n  "  '^o  "  '^  ^""^  '^ll  "  ^^^1'  "^12  "  ""l  ^  "^2'  ^^22  "  ^""l' 

6.  The  geometric  and  the  arithmatic  means  of  the  molecular 

radii  are  the  same, 

7.  The  geometric  and  harmonic  means  of  the  ionization  poten- 
tial are  the  same. 

8.  London's  formula  for  dispersion  forces  is  valid. 

M 
Thus,  the  simplified  expression  for  AE  can  be  written  as 


00 


r/d 


m 


Figure  1.   Radial  Distribution  Function, 
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AE^  =  27lSf^(N^V^  +  l^^V^)il/^T}/^ 
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L      1        .r^ 


^^ 


^2  ^    8r^ 


■/ 


1   ^  ■  ~fe  ;s(y)dy 


y         ny 


(94) 


or 


AE^  =   (N.J^   +  N^V^) 


-,V        1/2 
"l\ 


AS 


V 


1/2      2 


J  ^1^2 


(95) 


\^aere  AE,    is   given  by 


V 

AE^   =-E^   = 


~2 

2TrN  k 


3      J      V,*       ..."-Z   J 


8V  r"      ^ 
1    1 


g(y)<iy 


(96) 


Defining  a  solubility  parameter,  d,    by 


V  V 


(97) 


Equation  (95)  becoi:ies 


AE^  =  (N^^V^  +  ^^V^)[d^   -   5^]'^^^   ^^7//2 
The  partial  racial  energy  of  mixing  for  component  2  is  given  by 


(98) 


AE2  =  s^iV2C6^  -  &^]^  (99) 

To  derive  an  expression  for  AS^  Hildebrand  and  Scott  (22) 
made  the  following  assumptions. 

1.   The  internal  degrees  of  freedom  of  the  molecules  (rotation, 
vibration,  etc.)  are  the  same  in  the  mixture  and  in  the  pure  liquids 
as  in  the  dilute  gas. 
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2.    The  sane  free  volume  is  accessible  to  both  species  in 
the  mixture  and  this  is  equal  to  the  sum  of  the  free  volumes  of  each 
component . 

'•    !i    !i 

With  these  assumptions  the  Flory-Huggins  expression  (23)  for  partial 
molal  entropy  of  mixing  has  been  obtained. 

V, 


AS^  =  -R 


Ivep^+f^     (^^  "  V^)  J  (^°°^ 


If  V  is  assumed  to  be  equal  to  V  ,  we  get  the  ideal  entropy  of  mixing 

ASg  ~  -Rlnx^  (101) 

Using  equations  (100)  and  (101)  for  AS   in  equation  (9I)  the  following 
expressions  for  partial  molal  free  energy  of  solution,  AG„  are  obtained: 

—       2        2  ^2 

For  V^AV^,     AG^  =  V^-^^  (5^  -  62)   +  RTdn^^  +l^i^l  "  v   ^     ^'°^^ 

For  V^  =  V,,   AG2  =  V2^3^^(6^  -  6^)^  +  RTlnx^  (103) 


If  the  pure  component  2  at  the  pressure  and  temperature  of  the  solution 
is  chosen  as  standard  state, 


AG^  =  RTlnry^x^ 

where  7^  is  activity  coefficient  for  component  2  in  the  solution. 

Thus,  the  final  equations  for  a  2-component  solution  are 

V^     V^^(6,  -  6  )^ 
Irty^x^  =  Inip^  +  ijJ-^a   -  ^  )  +      ^: ^-  (104) 
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In-z^x,  =  Inx^  +-= — =—  (105) 

Equarions  (104)  and  (105)  have  been  used  extensively  for 
correlating  the  solubility  of  nonpolar  gases  in  nonpolar  solvents  by 
many  workers.   Among  those  are  Gjaldbaek  and  Hildebrand  (24,25)  who 
modified  equation  (104)  to  give  for  gas  solubilities 

i    /     "^2  \      ^2   ^^  2 

-Inx^  =  -Inx^  +  (  1  -  :^  W  In  :^  +  ^  (0^  -  6^r  (106) 

where  x^  is  ideal  solubility,  given  by  f„  =  f^x„ 
f^  =  fugacity  of  component  2  in  solution 

f_  =  fugacity  of  component  2  in  standard  state,  i.e.  pure  component 
2. 
In  obtaining  the  above  equation,  these  authors  used  the  condition 

N„  «  N  ,  and  thereby  set  ^  and  -if/     of  equation  (104)  approximately 
V^  -      1 

equal  to  x^  rp  and  1,  respectively.   Gjaldbaek  and  Anderson  (26),  and 

Gjaldbaek  and  Niemann  (27)  studied  the  solubility  of  nonpolar  gases  in 

polar  solvents.   They  included  the  empirical  factor  CO^  of  Hildebrand 

and  Scott  (22)  to  account  for  the  dipole  forces  of  the  solvent  and 

obtained  the  following  modified  expression: 

i   ^2          2     "^2   2 
^       I 2 .^ !_ 2_j  , RT_lj  ^^Q^^ 

A  B 

Gjaldbaek  and  Anderson  (26)  calculated  parts  (A)  and  (B)  of  equation 
(107),  and  found  that  the  difference  between  part  (A)  and  the  experi- 
mental values  of  lnx„  was  approximately  proportional  to  the  dipole 
m.oment  of  the  solvent.   However,  this  difference  was  not  equal  to  the 
part  (B) .   Thus,  they  found  that  equation  (107)  failed  to  predict  the 
solubility  of  nonpolar  gases  in  polar  solvents  correctly;  they  attributed 
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this  failure  to  the  unjustifiable  assumptions  made  in  the  derivation 
of  equation  (107).   Gjaldbaek  and  Niemann  (27),  and  Lannung  and 
Gjaldbaek  (28)  also  correlated  the  solubility  of  methane,  nitrogen, 
argon  and  ethane  in  alcohols  and  water  with  dielectric  constants  of 
the  solvent  and  the  polarizability  of  the  solute  gas.   However,  they 
did  not  find  good  agreement  with  the  experimental  values. 

Prausnitz  (29,30),  in  a  series  of  papers  on  gas-liquid  solutions, 
concer.aeu  that  since  regular  solution  equations  do  not  take  into 
account  the  large  volume  change  accompanying  the  dissolution  of  a  gas, 
they  cannot  be  used  without  modification  for  such  systems.   The 
required  modification  was  achieved  in  a  two  step  solubility  process. 
In  the  first  step,  the  gas  is  condensed  to  the  partial  molar  volume 
which  it  possesses  as  a  solute  in  the  liquid  solvent.   The  free  energy 
change,  AG^  for  this  process  is 

fj- 
AgJ  =  RTln  -^  (108) 

f 
^2 

where  f^  is  the  fugacity  of  the  hypothetical  liquid  at  one  atmosphere 
pressure. 

fj    is  the  fugacity  of  the  pure  gas  at  the  initial  condition, 
taken  as  one  atmosphere 

In  the  second  step,  the  hypothetical  liquid- like  fluid  was  dissolved  in 

the  liquid  solvent.   The  free  energy  change  AG2   for  this  step  is 

Ag"  =  RTln^/^x,  (109) 

—     —I    —II 
Using  the  condition  of  equilibrium  i.e.,  AG-  =  AG  +  AG^   =  0  for  the 

dissolution  process  and  the  regular  solution  equation  (105)  Prausnitz 
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obtained  the   following  relation  for  solubility  of  a  nonpolar  gas  in 
nonpolar  solvents: 

L        L      ■    \    - 


2 

tor  tne  hyp 


To  deter^iine  the  parameter  f^  for  the  hypothetical  liquid  state,  he 

.1 

hypothesized  thac  the  reduced  fugacity,  _2  ,  is  a  universal  function  of 

p 
'c 

r-,  -r 

reduced  tem-perature  ~    •      To  obuain  V^ ,  it  was  assuined  that  its  value 

c 

is  independent  of  the  solvent.   The  author  simultaneously  estimated  the 
3  parameters  fi",  V"  and  6^  from  the  ey.isting  solubility  data  for 
the  particular  gas  in  various  advents  by  fitting  rhe  data  to  equation 
(110).   The  fitting,  lo  this  equation  was  done  in  such  a  way  that  the 
values  of  fT  ooeyed  the  corresponding  scaces  lav;  when  plotted  as 
z'^l^ ^   against  reduced  temperature.   lie  found  ;;hat  equation  (110) 
correlates  the  available  solubility  daca  quite  well  for  most  gases 
except  hydrogen. 

Prausnitz  and  Shair  (31),  in  an  a'tcoimpt  to   correlate  zhe 
solubility  of  nonpolar  gases  in  polar  solvents,  rewrote  equation  (110) 
in  the  following  form: 


^7,  (111) 


Their  hypochesis  was  that 

L 
lir/,^  =  V^  ?(T,6.  ,  properties  of  solvent)  (112) 

and  the  parameters  V.T*,  5^_,  and  f]!^  c^t^-j aiA   only  upon  solute  properties; 
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therefore,  these  could  be  determined  from  the  solubility  of  the  given 
solute  gas  in  other  solvents.   Thus,  they  could  correlate  gas  solubility 
data  in  polar  solvents  by  plotting  isotherms  of  (ln7^)/vr'  vs.  6„ . 

Yen  and  McKetta  (32)  have  presented  a  correlation  of  nonpolar 

gas  solubilities  in  polar  nonassociated  liquids.   They  used  the 

M 
following  basic  equation  for  AE  ,  which  was  obtained  by  setting 

g..  =  g.  in  equation  (92): 


—  /  O  or  dr 
„2  /  ^22°2 
V  ^ 
2  0 


(113) 


However,  the  authors  considered  the  electrostatic  interaction  among 
the  polar  molecules  by  using  the  Stockmayer  potential,  i.e.. 


'i^''    /fl^^.ill 


2 

(114) 


^i=^^i[(t)   -  [^  )     -7^(^.^) 


where  'i>^^    is  potential  energy  of  interaction  of  polar  molecules 
e.    =  Lennard- Jones  minimum  potential  energy 
T}(6,0)  =  angular  dependence  of  dipole-dipole  interaction 

It  was  shou^  that  induction  forces  between  the  nonpolar  and  the  polar 
molecules  were  small,  and  these  were  neglected.   They  made  the  same 
basic  assumptions  as  those  used  by  Hildebrand  and  Scott  (22)  in 
obtaining  equation  (105),  and  obtained  the  following  equation  for  the 
activity  coefficient  of  the  solute  gas  in  the  solution: 


wnere 


1^1 
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and  is  a  contribution  of  the  dipole  interaction  of  solvent  molecules. 

For  the  solubility  of  a  nonpolar  gas  in  a  nonpolar  solvent, 
A  =  0  and  equation  (115)  reduces  to  equation  (105). 

Yen  and  McKetta  used  the  two  step  process  of  Prausnitz  and 
Shair  (31),  described  earlier,  to  account  for  the  volume  change  on 
mixing.   A  procedure  similar  to  that  given  in  reference  (31)  v/as  used 
to  obtain  the  parameters  V_ ,  6^  and  f   for  the  hypothetical  liquid 
state.   It  was  shown  that  the  values  of  these  parameters,  as  estimated 
from  the  solubilities  of  nonpolar  gases  in  nonpolar  solvents,  could  be 
used  to  predict  the  solubility  of  nonpolar  gases  in  polar  solvents. 
Thus,  they  obtained  the  factor  A  as  a  characterizing  constant  for 
each  polar,  nonassociated  liquid.   It  was  found  that  the  value  of  A 
increases  with  increasing  solubility  parameter  of  the  liquid. 

In  a  recent  paper  Loeffler  and  McKetta  (33)  presented  a  thermo- 
dynamic correlation  of  nonpolar  gas  solubilities  in  alcohols  as  an 
extension  of  the  correlation  developed  by  Yen  and  McKetta. 
2.2.3    Cavity  Models: 

Cavity  models  have  been  proposed  in  an  attempt  to  explain  the 
solubility  of  gases  in  liquids.   Since  most  of  the  earlier  cavity 
models  are  approximate  in  nature,  they  are  only  briefly  mentioned. 
The  model  of  Reiss  _et  al.  (34),  based  on  the  ideas  of  the  scaled 
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particle  theory  (35)  and  extended  by  Pierotti  (36,37)  appears  to  be 
more  exact  and  has  therefore  been  considered  in  more  detail  in 
Section  2.2.4. 

It  v;as  Uhlig  (38)  who  first  considered  the  solubility  process 
to  be  made  up  of  two  steps:  (1)  creation  of  a  spherical  cavity,  for 
which  work  is  done  against  the  surface  tension;  (2)  introduction  of  the 
solute  molecule  into  the  cavity,  which  involves  energy  terms  arising 
out  of  interactions  between  solvent  and  solute  molecules.   Considering 
these  energies  he  obtained  the  following  expression: 

where  r]  =   Ostwald  solubility  coefficient 

a  =   surface  tension 

E  =  energy  of  interaction 

r  =  radius  of  tne  solute  molecule. 
El ley  (39,40,41)  accounted  for  the  energy  and  entropy  of  solution  of  an 
inert  gas  molecule  in  water  and  organic  solvents  by  using  the  same  two 
consecutive  steps  as  considered  by  Uhlig  (38);  however,  he  calculated 
the  energy  and  entropy  of  formation  of  the  cavity  by  the  following 
approximate  thermodynamic  relations: 

AS   =(a//3)V 
AE^  =Cra//5)v 


where   a  =  thermal  expansion  coefficient 
/3  =  compressibility  coefficient 


Karsaz  and  Halsey  (42),  and  Volk  and  Halsey  (43)  also  assumed 
the  solution  process  to  be  made  up  of  the  two  following  steps:   A 
vacancy  appropriate  to  contain  the  solute  molecule  is  first  formed  in 
the  solvent;  che  partial  molal  Gibbs  free  energy  of  making  this  void 
is  G  =  H  -  TS  .   Secondly,  the  void  is  filled  with  a  solute  molecule 

V      V       V 

and  allowed  to  mix  with  the  solvent.   They  considered  the  Lennard- 
Jones-Devonshire  model  in  which  the  solute  is  in  the  center  of  a  cell, 
surrounded  by  solvent  molecules.   3y  considering  the  free  energy  of  the 
system  in  the  gas  and  liquid  phases,  these  authors  obtained  the  follow- 
ing expression: 

P9 
log  -=^  =  CO  /kT  +  G  /kT  -  log(v/kT)  (116) 

"X^O         V  ov/ 

where   G)  =  minimum  potential  energy 

V  =  free  volume. 

Equation  (116)  cannot  be  immediately  used  because  the  theory  does  not 

give  directly  either  G  ,  CO  ,  or  v. 

V    o 

Kobatake  and  Alder  (44)  proposed  a  Lennard- Jones- Devonshire 

form  of  the  gas  cell  potential  similar  to  that  of  Karsaz  and  Halsey. 

2.2.4    Scaled  Particle  Theory  of  Reiss,  Frisch,  Eelfand  and  Lebowitz 
(34,35) 

This  theory  is  a  refined  version  of  the  cavity  models  discussed 
in  the  previous  section.   However,  it  is  discussed  separately  here 
since  it  forms  the  basis  of  the  theoretical  approach  used  in  the 
present  work.   The  principal  advantage  of  scaled  particle  over  previous 
theories  of  this  type  is  in  calculating  the  free  energy  of  cavity 
formation  in  a  reasonably  rigorous  fashion. 

Consider  a  mixture  consisting  of  components  1,2, m  (the 
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solvent).   Into  this  mixture  a  soluue  molecule  a  is  introduced.   The 

r^al  molecules  are  assumed  to  possess  hard  cores  of  diameter  a  

^m'^a'      ■^^"^^i"'g  that  the  total  potential  energy  is  the  sum  of  pair 
potentials,  the  following  expression  for  chemical  potential  of  the 
solute  can  be  obtained  (45). 

m     !°    ^ou   (r,4) 

,     1/2 
\   27imkT  ; 


where   A  = 


p.  =  number  density  of  component  j 

u   (r,^)  =  pair  wise  potential  between  solute  and  solvent  component  j 

g  . (r,^)  =  radial  distribution  function  between  solute  and  solvent 
component  j 

E,   =  a.   charging  param.eter  (see  ref .  45)  which  allows  the  solute 
molecule  to  be  coupled  with  the  solvent. 

We  note  that  only  interactions  between  the  solute  a  and  solvent  of 

species  j  are  involved  in  this  equation.   In  equation  (117)  the  lower 

and  upper  integration  limits  of  ^  correspond  to  complete  uncoupling 

and  complete  coupling  of  the  solute  molecule,  respectively. 

The  pair  potential  u  .(r,^)  is  assumed  to  be  of  the  following 

form  (see  Figures  2  and  3) 

u  .(r,«)  =  u^.  (r,t)   +  u^.(r,4  )  (118) 

aj       aj   ^h    aj  ^s 


For  convenience,  we  use  subscript  a  to  denote  solute  m.olecule 
and  subscripts  1,  ....  m  to  denote  solvent  species.   This  nomenclature 
is  not  the  same  as  used  in  the  previous  sections. 
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Hard  Core 


"C^-Pf" 


Sort"  Outer  Part 


Figure  2.   Schematic  Diagram  of  Molecule  of 
Species  j . 
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'igure  3-   Pair  Potential  Between  Solute  and 
Solvent  Molecule  of  Species  j. 
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where  u  .(r,£,  )  and  u  .(r,e  )  are  respectively  the  hard  sphere  and 
aj   ^n      aj   ^s 

soft  portions  of  the  potential  between  solute  and  solvent  component  j . 
^  and  ^  are  charging  parameters  for  the  hard  and  soft  potentials, 

respectively . 

a  +  a . 
u''.(r)=co         r<^a.=-^T ^  (119) 


<J^^^=° 


u^(r,^  )  =  i   u  .(r) 
CL2         s         ^s  aj 


r  >  a  . 


By  allowing  ^,  to  vary  from  0  to  a  .  we  "charge  up"  the  hard 

part  of  the  potential.   The  term  |  is  the  'scaling'  parameter  (hence 

the  name  'scaled  particle  theory').   Thus  when  ^,  =  a^  .  the  hard  core 

of  solute  a  is  fully  coupled  to  the  rest  of  the  system.   ^  varies  from 

0  to  1.   When  |  =  1,  the  soft  part  of  the  potential  will  be  fully 
s 

coupled.   A  value  of  Zc;ro  assigned  to  both  ^,  and  ^   implies  that 
'^  h      s 

solute  a  is  decoupled  from  the  system. 

The  solute  chemical  potential  may  now  be  written 

3    ^    s 
u  =  kTlnp  A-^  +  g  -rg^  (120) 


where 


h    ™    p      no'^r.^^^'^'J 


/■      r.       r=0 
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g.  =  -P.  /  -^y.     I  ^ ^aj^^'^h'^s  =  ^^"^^^^^ 


(121) 


^h=0 

s   m   ^Sp 
i^  =  So.  /  d^   /  u^(r)g..(r,e,_  =  a_,e.)47lT"dr       (122) 


j"^    e'=o    a'  . 


^s     aj 
The  charging  procedure  involved  in  equation  (121)  is  tantamount 
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to  first  introducing;  a  hard-sphere  solute  molecule  of  diameter  a  . 

•=  a 

In  evaluating  g"  a  solute  molecule  of  diameter  (-a^)  is  introduced, 
°  ^a  i^ 

and  then  scaled  up  to  its  full  diameter  a  . 

For  a  vapor  and  liquid  phase  in  equilibrium 

u.^   =  u^  (123) 

^a.      •  a. 

G      L 
where  l;  and  u  are  chemical  potentials  of  the  solute  molecule  in  the 

G     , 
gas  and  liquid  phases,  respectively.   The  chemical  potential,  u^  ,can  be 

written  (45): 

3 

Q 

where  f   is  the  fugacity  of  the  solute  m.olecule  in  the  gas  phase, 
a 

By  equating  (12C)  and  (124) 

^G      _h   _s 

The  mole  fraction  is 

a   ^  p. 


.G      _h   _s 


(126) 


or   in  terms   of   partial  m.olal   quantities 

n  _h       _s 

a  \  a    ,     a    ,    ,    ,  RT 


a  /  J 


where   2  n      =  number   of  mioles    in  volume  V. 
From  equation   (10),    we  have 


ka 


f 

—  =  7  K 

X    'a  a 
a. 

where  7  is  the  activity  coefficient  of  solute  a  in  the  solution  and 

K  is  the  Henry's  law  constant  defined  as: 


lim  f   =  K  X 


X  -*■  0 
a 


X.  =  0 
J 

J  *   1 

Thus  equation  (126)  can  now  be  written  as: 

_h   _s 

^a  ^  ^a  ^  ,  ,   RT 


l^Va>=l^-^k¥^^^(-^2"j)  (128) 

J 


Evaluation  of  g  : 

^ 

By  applying  the  treatment  of  Reiss  et_   al .  (35)  to  a  mixture  containing 
m  solvent  components  and  assuming  that  the  concentration  of  solute 
is  very  low  so  that  u   =  0.  it  can  be  shown  that 

g   =  2  p   /  G   (r)47n:  dr  (129) 

"■  j=#a  J  '   "-J 
0 

where  G  .(r)  is  the  radial  distribution  function  for  solvent  component 
j  around  a  solute  molecule  a  when  the  two  are  in  contact.   The  quantity 

on  the  right  hand  side  or  equation  (129)  can  be  shown  to  be  equal  to 

a 
a 
the  work  of  introducing  a  cavity  of  radius  — r  into  the  solvent,  which 
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we  term  W(  —~   ).   The  word  cavity  here  is  used  to  mean  a  spherical 

volume  from  which  all  hard  core  parts  of  solvent  particles  of  all 

species  are  excluded;  i.e.,  the  cavity  does  not  imply  only  exclusion 

of  solvent  particle  centers. 

a 
Lebowitz  et  al.  (46)  have  obtained  an  expression  for  W(  —  ) 

from  -probability  considerations,  as  outlined  below. 

From  the  relation  between  the  probability  of  a  configuration 

and  the  work  of  creating  that  configuration,  we  have  (47) 

a  a 

W(  ^  )  =  -kTlnP^(  -^  )  (130) 

a 
where  Pn(  ~^  )  is  the  probability  that  there  are  no  parts  of  any 

a         a  a 

solvent  particle  in  the  cavity  of  radius  -^  .   For  — j  <^  0,  W(  -^  )  is 

a 
obtained  as  follows.   Consider,  for  example  the  case  of  -^  =  0,  i.e., 

a  point  solute  particle  located  at  a  point  X.   The  nearest  the  center 

of  a  j  solvent  particle  can  come  to  X  is  a./2.   The  probability  that 

particle  1  of  species  j  (a  particular  particle)  is  in  a  spherical 

volume  of  radius  r  is: 

4   3 

3^ 


^  3 

47nr  N. 

and  the  probability  thac  any  j  particle  is  in  this  space  is  ^ 


3V 

For  a  sphere  of  radius  a./2  this  becomes 

J 

3  7r(  ^  )  p^ 

Similar  probabilities  apply  for  other  species.   The  total  probability 
that  a  molecule  of  any  type  is  in  the  sphere  of  radius  a./2  is: 


It  should  be  noted  that  the  shell  radius  a./2  referred  to  varies  with 
the  species  involved. 

The  probability  that  a  cavity  of  radius  zero  exists  is  then 


Pq(0)  =  l.fTr2  (/;P.  (131) 

a 
The  same  arguments  apply  for  any  value  of  — ^J^  0,  so  that 

,m/a.\3       a 
PqCI  )  =  1  -f  rSJ^-^  )   p.    ^<  0        (132) 

using  equation  (130), 


m   /a  +  a  \  3 

i 


2      )  P.  ■ 


J  - 


a 

-f  4:  0         (133) 


From  thermodynamic  considerations,  for  very  large  r,  the  following 
expression  for  the  work  of  formation  of  a  cavity  of  radius  r  can  be 
written  (34,35) 


W(r)  =  -  Trr  p  -r  Attt  a^   1 

J  0  I      r 


+f(p.,T)  (134) 


where  P  is  the  pressure,  cr  is  the  interfacial  tension  between  the 
fluid  and  a  perfect  rigid  wall,  and  6  is  a  distance  of  the  order  of 

the  thickness  of  the  inhomogeneous  layer  near  the  interface  and  accounts 
for  the  curvature  dependence  of  the  surface  work.   f(p.,T)  is  an 
arbitrary  function  of  density  and  temperature,  not  dependent  on  r. 
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For  r  positive  but  small  W(r)  is  expanded  in  a  Taylor  series  about  r  =  0, 


W(r)  =  W(0)  +  rW  (0)  -f  ^  ^~W"(0)  +  K  ^^ 


r  >  0 


(135) 


Since  W(r)  and  its  first  two  derivatives  have  been  shown  (34) 
to  be  continuous  at  r  =  0,  the  first  three  coefficients  of  the  series 

can  be  evaluated  by  coiriparison  with  equation  (133).   The  coefficient 

3 
of  the  term  in  r   can  be  obtained  by  comparison  with  equation  (134). 

Thus 


W(^) 

kT 


r^^: 


■1^(1  -  C3)  -^  ,  ,  . 


a 


r  12^. 


1  -  C. 


ist? 


(1  -  C3) 


(2)     ^I^fe(-f 


where 


m 


Ci  =t7r2p(a  )■ 
2=1 


(136) 


— s 

Evaluation  of  g  ; 


m 


e  =1 


J^   ^  =0   a  . 


a  .  ,i;  )47ir  dr 
aj   s 


s 

g   can  be  decomposed  as  follows: 


— s   _s   _^s    _s 

g   =  e  -  is   +  pv 
a   a    a   ^  a 


(137) 


g   2 

where  e^ ,  v  and  s,  are,  respectively,  internal  energy,  volume  and 


entropy  associated  with  the  charging  of  the  soft  potential, 
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s 

At  atmospheric  pressure  the  terni  in  pv   can  be  shown  to  be 

s 

small  and  can  be  neglected  in  comparison  to  "e   (36,37).   The  entropy 

,  g 
of  charging  s   is  negative  because  the  system  becomes  more  organized 

upon  charging.   Since  it  is  difficult  to  make  a  quantitative  estimate 

of  this  quantity,  it  will  be  neglected  as  a  first  approximation. 

g 
e  can  be  obtained  as  follows:   Consider  a  shell  of  thickness  dr  at  a 

a 

distance  r  from  the  solute  molecule. 

The  number  of  particles  of  solvent  component  j  in  this  shell  is 

given  by 

2 
(47:r  dr)p  g  (r) 
J  '^J 

2   can  now  be  ^^7ritten  as 

00 

i^i^e  =  .2  /   u  .Ato  p.g  .(r)dr  (138) 

a  . 

aj 

Evaluation  of  g  by  the  above  equation  requires  a  knowledge  of 

cx 

g  .(r),  which  is  difficult  to  obtain  for  the  systems  studied  here. 

As  a  first  approximation  it  will  be  assumed  that  solvent  particles  are 

uniformly  distributed  around  the  solute  molecule  for  r  larger  than 


a  .  I.e. , 


Thus  m 


g  .(r)  =1  r  >  a  . 


g%  2p.  \   V.^Tir^dr  (139) 


Interaction  Energy,  u  . 
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The  interaction  energy  of  the  nonpolar  solute  molecule  in  an 
aqueous  potassiuir.  hydroxide  solution  can  be  written  as  the  sua  of  the 
individual  pair  contributions  u  .,  u^ „  and  u  „  where  subscripts  1,  2 
and  3  denote  water,  K  and  OH  . 

If  contributions  due  to  quadruple  moments  are  neglected  these 
interactions  are: 

no 
where   u  *  =  nonpolar  interaction 

OL  — 

u  C"'   '^     =   induced  interaction  between  nonpolar  solute  molecule  and 
the  permanent  dipole  of  a  water  molecule. 

_  np  ,   (C.indu)       .    „  ^  /i/-,n 

u.=u.-u.    ■         1=  2,3  (141) 

ai    ai    ai  ' 

where  ji    .    =   induced  interaction  between  the  nonpolar  solute  molecule 

and  the  electrical  charges  on  ion  i. 

The  nonpolar  interaction  energy  between  the  solute  molecule 
and  the  solvent  components  is  made  up  of  dispersion  and  repulsive 
energies.   It  V7ill  be  assumed  that  for  the  systems  studied  here  this 
contribution  is  adequately  represented  by  the  Lennard- Jones  6-12  poten- 
tial. 


r/J.\12  .a.\6 

np        .  if     ai    \  /     aj 

ctj     aj  [^  V  r         \  '^ 


(142) 


The  follov7ing  relation  is  assumed  between  pure  component  Lennard- Jones 
parameters  and  those  in  a  mixture. 

1/2  ^c   ^  ^i 

£  .  =  (£  £.)^^  C    .     =-^^ 1  (143) 

The  interaction  between  a  single  permanent  dipole  and  the 
induced  dipole  in  a  neutral  solute  molecule  after  statistical  averaging 


^h 


over  all  possible  orientations  is  (4S): 

9 


al     = T^  (144) 

r 

where  u.,  is  the  permanent  dipole  moment  of  water  molecule  and  p  is 

1  a 

the  polarizability  of  the  solute  molecule.   In  the  present  case,  a 
solute  molecule  is  surrounded  by  not  one  but  several  water  dipoles; 
therefore  at  any  instant  of  time  the  dipole  induced  due  to  all  these 
dipoles  should  be  considered.   The  exact  calculation  of  this  contribu- 
tion can  be  made  only  if  the  distribution  of  the  water  molecules  around 
the  neutral  solute  molecule  is  known.   However,  in  this  model  the 
dipoles  are  assumed  to  be  freely  rotating  so  that  we  can  average  over 
all  orientations  of  each  dipole.   Consequently,  equation  (144)  for  the 
dipole- induced  dipole  interaction  can  be  used  even  in  the  present 
situation. 

The  following  assumptions  are  implicit  in  equation  (144)  for 
the  induced  interaction  (49). 

1.  The  solute  and  water  molecules  are  electrically  neutral. 
The  permanent  dipole  of  the  wacer  molecule  is  characterized  by  two 
equal  charges  of  opposite  sign. 

2.  The  distance  r  between  molecules  is  large  compared  to  the 
separation  of  charges  within  each  molecule. 

3.  The  field  produced  by  each  molecule  at  the  other  molecule 
is  homogeneous . 

4.  The  solute  molecule  is  polarizable  and  possesses  isotropic 
polarizability. 
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5.  Each  dipole  moment  is  that  in  zero  external  field. 

6.  The  effective  mutual  potential  energy  of  the  pair  of 
molecules  is  the  ensemble-average  over  all  relative  orientations. 

The  ion-ir.duced  dipole  interaction  is  given  by 

(Ci.indM)  ^  _   r  (ind)^^ 
0 

where  c  is  the  electric  field  produced  by  the  ions  considered  at  the 

-■,         -,-.-,^.1,       r-    •         ■         1      (ind) 
position  or  the  neutral  molecule,  C.  is  the  charge  or  ion  i  and  u 

is  the  dipole  induced  in  the  neutral  solute  molecule  due  to  the  ions. 

-1  -  ■  -,  J     ^      n    (ind)  .    .    , 
Assuming  external  riexds  to  oe  small,  /i      is  given  by: 

In  the  case  of  a  single  ion 

e . 

r 

(e .  ,  xndji)  p  e . 

U."^       =2/p— rdr 
^ai  /  ^a  :> 

^'     r 

~   2 

=  -  ^T^  (147) 

2r" 

In  the  present  situation,  the  field,  ^  is  not  due  to  a  single 

ion,  but  to  all  the  ions  present  and  therefore  its  correct  value  can  be 

obtair.id  only  if  the  distribution  of  the  ions  around  the  neutral 

molecule  is  known.   However,  we  have  assumed  that  the  ions  are  uniformly 

distributed  around  the  neutral  molecule.   In  other  words  there  is,  on 

the  average,  a  spherically  symjnetrical  charge  distribution  around  the 
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neutral  molecule.   For  such  a  charge  distribution,  the  electric  field 
at  any  point  in  the  neutral  molecule  can  be  obtained  as  follows  (50) . 

Consider  a  point  P  in  the  neutral  molecule  and  imagine  a 
thin  cone  of  solid  angle  cCj   and  vertex  at  P  to  cut  the  sphere  in 

areas  oS^    and  6S   (Figure  k) .   The  right  section  of  the  cone  at  6S 

2 
is  r'TiSc  and 

2  2 

r  do)  r..6co 

T-:; —  =  cosa,  or  63^  =  

Oa,  1   cos 


Similarly  p 


6S^  = 


rp5:o 


2   cosa 


and  if  cr  is  the  surface  density  of  charge  upon  the  sphere,  the  charges 
upon  6S..  and  5S   are 

2  2 

Cr,6a)       cr^Sco 

and 


cosa        cosa 


The  electric  field  at  point  P  is  given  by 


2c      2. 

r,5a)   or  Oco 

'       -  -  =  0  (li^) 


2       2 
r  cosa  XpCosa 


Thus  v;e  see  that  for  a  spherically  symmetrical  charge  dis- 
tribution around  the  neutral  solute  molecule,  the  electric  field 
produced  at  any  point  in  the  solute  molecule  is  zero.   Consequently 
the  induced  dipole  introduced  by  all  the  ions  acting  together  is 
zero;  hence  there  is  no  contribution  due  to  the  ion-induced  dipole 
interaction. 

Although  the  above  argument  is  valid  for  such  an  average 
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Figure  k.      A   Sphere  Having  a  Uniform 
Surface  Density  of  Charge. 
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situation,  at  any  instant  we  should  expect  for  a  real  solution  that 
a  small  finite  fluctuating  field  £"  exists  at  the  point  P  in  the 
solute  molecule,  so  that  at  such  an  instant  the  ion- induced  dipole 


interaction  is 


(C.,ind,Li) 

U    . 
'^CXl 


f  'c£  '^  ■ 


ii) 


(149) 


^0 


This  contribution  to  the  intermolecular  potential  will  always  be 


negative  (an  attractive  force)  for  ^'  /  0  regardless  of  the  sign 
of  ^  .      The  contribution  indicated  by  equation  (149)  is  presumably 
small,  and  is  neglected  in  this  work. 

Thus  the  potential  energy  of  interaction  of  the  solute  molecule 
is 


.J      r/u.\o     /u.\ 
2u.=2  4e.    (-^)  -     ( -^] 


a  .  xl2  1    P^Mi 


(150) 


Substituting   equation   (150)    in  equation   (139),    we  obtain 


— s 
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..r 


=  -1671  /  2  p.e   . 


ai 


j=l  J  ^J   l-c 
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=    -loTT         2  p.£     . 
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At   a 


aj 


C    . ,    e     can  be  written  as: 
aj       a 


_s 

e      = 
a 


32^  V       ^3 
9    ,^,^j  a:  aj 


^1    g.l 

al 


(152) 


CHAPTER  3 
EXPERIMENTAL  APPARATUS  AND  PROCEDURE 

A  variety  of  approaches  have  been  used  to  determine  the 
solubility  of  gases  in  liquids.   A  recent  review  on  the  solubility 
of  gases  in  liquids  by  Battino  and  Clever  (4)  describes  these  methods 
in  detail.   Consideration  of  the  existing  methods  of  measurement  of 
gas  solubility  reveals  that  while  there  are  many  approaches  to 
measuring  solubilities  precisely,  only  a  few  methods  are  suitable  for 
determining  very  small  amounts  of  dissolved  gases.   In  particular  the 
very  accurate  manometric  methods  (56)  cannot  be  applied  when  the 
solubility  is  as  low  as  those  experienced  in  the  present  work.   Gas 
chromatography  is  particularly  suitable  for  the  measurement  of  low 
solubilities;  moreover,  it  is  rapid  and  has  the  necessary  versatility. 
Sample  size  can  be  varied  to  a  large  extent,  and  the  dissolved  gas 
can  be  stripped  from  the  solution  and  concentrated  by  a  suitable 
technique  before  analysis. 

Paglis  (51)  and  Elsey  (52)  have  employed  this  technique  for 
the  determination  of  dissolved  oxygen  in  liquid  petroleum  and  lubricating 
oils.   Swinnerton  (53)  described  a  method  for  the  determination  of 
dissolved  gases  in  aqueous  solutions  based  on  gas  chromatography. 
Gubbins,  Garden  and  Walker  (54)  m.ade  important  modifications  in  the 
method  of  Swinnerton  (53)  and  obtained  gas  solubilities  in  electrolytic 
solutions.   Most  of  the  data  reported  here  were  obtained  by  the  method 
of  Gubbins,  Garden  and  Walker  (54)  (henceforth  called  Method  I), 
whereas  some  of  the  data  for  systems  exhibiting  very  low  solubilities 
were  obtained  by  a  modification  of  this  method  employing  a  concentration 
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technique  (Methods  II  and  III) . 
Method  I 

A  Perk in -Elmer  gas  chromatograph  equipped  with  thermal  conducti- 
vity and  flame  ionization  detectors  was  used  in  this  work.   The  output 
signal  from  the  detector  was  recorded  on  a  1  m.v.  Sargent  recorder 
(Model  SR)  provided  with  a  disc  integrator  (Model  204).   All  analyses 
were  made  using  the  thermal  conductivity  detector.   Helium  was  used  as 
the  carrier  gas  for  all  solute  gases  except  hydrogen  and  helium;  for 
these  nitrogen  was  used  as  the  carrier  gas.  A  gas  flow  rate  of  approxi- 
mately 30  ml.   per  minute  was  used  in  all  the  measurements. 

The  solution  sample  to  be  analyzed  was  injected  into  a  glass 
stripping  cell  outside  the  chromatograph  (Figure  5).   The  carrier  gas 
flowing  through  the  sample  line  of  the  chromatograph  was  diverted 
through  the  stripping  cell  and  its  associated  equipment  before  re-enter- 
ing the  sample  column.   The  glass  cell  had  a  medium  porosity  glass 
frit  (fine  frits  were  found  to  give  rise  to  excessive  pressure  drop 
when  wetted)  and  was  equipped  with  a  rubber  septum  for  the  injection 
of  the  sample,  and  a  teflon  stopcock  to  facilitate  draining  of  the  cell. 
Gubbins,  Garden  and  Walker  (54)  found  that  a  cell  having  a  volume  of 
7  ml. was  most  satisfactory  for  general  purposes.   In  this  work,  however, 
attempts  were  made  to  increase  the  sensitivity  of  the  method  by 
increasing  the  size  of  the  stripping  cell  to  accomodate  larger  samples. 
Increasing  the  diameter  of  the  cell  led  to  excessive  broadening  and 
tailing  of  the  chromatograph  peaks,  apparently  because  of  mixing  in 
the  gas  space  above  the  liquid.   It  was  found  that  a  stripping  cell 
about  1  cm.  in  diameter  and  having  a  height  of  70  cm.  could  be 
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successfully  employed  to  increase  the  sensitivity  of  the  method.  It 
v;as  later  found  that  Swinnerton  (55)  had  also  used  a  similar  type  of 
glass  chamber  in  his  work. 

A  solution  sample  of  known  volume  (20-40  ml.)  containing  the 
dissolved  gas  was  injected  through  the  rubber  septum.   The  carrier 
gas  was  dispersed  by  the  fritted  glass  disc  into  a  stream  of  fine 
bubbles  and  stripped  the  dissolved  gases  from  the  solution  quickly  and 
completely.   Before  returning  the  gases  to  the  chromatographic  column 
they  were  passed  through  two  drying  columns  containing  Drierite  to 
eliminate  water  vapor,  followed  by  a  stainless  steel  coil  immersed  in 
a  water  bath  to  bring  the  gases  to  ambient  temperature  before  return- 
ing them  to  the  column.   This  was  particularly  important  when  deter- 
mining solubility  at  very  low  or  very  high  temperatures,  since  the 
sudden  cooling  or  heating  of  the  gas  stream  otherwise  caused  base- 
line instability.   Care  was  taken  to  keep  all  tubing  of  small  diameter 
to  prevent  back  mixing  and  broadening  of  the  peaks.   The  drying  columns 
were  1  en;,  in  diameter  and  30  cm.  in  length.   The  use  of  rubber  tubing 
v;as  completely  avoided  to  minimize  contamination  by  air  diffusing 
through  it. 

This  technique  of  introduction  of  the  gas  sample  by  stripping 
from  a  solution  results  in  slightly  wider  and  less  symmetrical 
chromatograph  peaks  than  those  obtained  when  a  gas  sample  is  injected 
directly  into  the  chromatograph  sampling  port.   However,  it  does  not 
seem  to  result  in  any  loss  of  accuracy  or  precision.   The  amount  of 
the  dissolved  gas  was  obtained  by  means  of  an  appropriate  calibration 
factor. 
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For  these  gases  whose  solubility  in  water  is  accurately  known 
at  25  C,  the  calibration  factor  was  obtained  by  carrying  out  the  same 
procedure  on  a  water  sample.   VJhen  ;his  information  was  not  available, 
an  accurately  known  volume  of  pure  dry  gas  was  injected  directly  into 
the  chromatograph  to  serve  as  a  standard.   The  relationship  of  the 
amount  of  dissolved  gas  to  the  instrument  response  was  checked  and 
was  found  to  be  linear  within  the  limits  of  the  experimental  error. 

Large  discrepancies  often  exist  between  the  solubility  values 
reported  by  different  workers.   Cook  and  Hanson  (56)  have  listed  the 
following  factors  which  need  careful  attention  in  order  to  minimize 
these  discrepancies: 

1.  Complete  degassing  of  the  solution 

2.  Ascertaining  the  true  amounts  of  gas  dissolved 

3.  Attaining  equilibrium 

4.  Making  certain  that  transfer  of  the  gas  from  the  primary 
container  to  the  apparatus  does  not  involve  contamination. 

Each  of  these  factors  will  be  discussed  with  particular  reference  to 
the  procedure  adopted  in  this  work. 
Degassing  the  Solution: 

Solutions  were  prepared  from  specially  distilled  (all  glass- 
teflon  still)  and  degassed  water.  The  method  of  degassing  water  was 
similar  to  that  employed  by  Clever  et  al.  (57),  and  consisted  of  two 
stages.  The  first  stage  involved  boiling  the  water  to  evaporate  a 
portion  of  it  and  to  remove  perhaps  907o  of  the  dissolved  gases.  Then 
in  the  second  stage  this  substantially  degassed  solvent  was  agitated 
with  a  magnetic  stirrer  under  vacuum.   Pouc^sium  hydroxide  pallets 
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used  for  preparing  aqueous  solutions  were  Baker  analyzed  reagent  grade 
and  contained  a  maximum  of  1%   K^CO^ .   Bottles  containing  potassium 
hydroxide  solution  were  fitted  with  glass  bulbs  containing  ascarite, 
so  that  air  entering  the  bottle  at  the  time  of  withdrawal  of  solutions 
was  free  from  carbon  dioxide.   Solutions  prepared  in  this  way  when 
injected  into  the  chromatograph  after  equilibrialion  with  the  carrier 
gas  gave  no  peak. 
Saturation  Method: 

Saturated  solutions  of  the  gas  in  the  electrolyte  were  prepared 
by  bubbling  the  gas  by  means  of  a  medium  porosity  glass  frit  through  a 
set  of  presaturators  and  then  through  the  saturating  vessel.   The 
saturating  vessel  was  a  3  neck,  500  ml.  round  bottom  flask.   The  central 
neck  accomodated  a  teflon  clad  ground  glass  fitting  in  which  in  turn 
was  installed  the  glass  frit  for  bubbling  the  gas.   The  other  two 
openings  were  used  to  withdraw  samples  and  as  an  exit  for  the  gas.   The 
gas  leaving  the  saturator  bubbled  through  a  glass  test  tube  containing 
the  same  solution  as  in  the  saturating  vessel  to  avoid  any  contamination 
due  to  atmospheric  air.   Samples  were  removed  by  means  of  a  syringe 
fitted  to  a  teflon  needle  dipping  in  the  electrolyte  by  means  of  a 
Hamilton  teflon  valve. 

Attainment  of  equilibrium  was  checked  by  withdrawing  samples  after 
different  intervals  of  saturation  time  and  injecting  into  the  stripping 
cell.   The  equilibrium  solubility  vjas  taken  as  the  value  measured  when 
at  least  3  samples  withdrawn  at  different  times  gave  the  same  peak  area. 
This  process  was  repeated  at  different  gas  flow  rates  and  it  was  found 
that  gas  flow  rate  did  not  affect  the  equilibrium  solubility  value. 
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However,  to  avoid  ioaming  in  the  presaturators  and  saturating  vessel, 
a  gas  flov7  rate  of  about  30  xl.  per  min.  was  used.   At  this  flow  rate  it 
took  about  25  min.  to  attain  equilibriuni. 

The  possibility  of  supersaturation  was  guarded  against  by  the 
following  procedure.   A  sample  of  the  solution  through  which  solute 
gas  was  bubbling  was  removed  and  analyzed.   Bubbling  of  the  gas  was  then 
stopped,  and  the  KOH  solution  maintained  under  an  atmosphere  of  the 
solute  gas.   Samples  of  the  solution  were  analyzed  periodically,  and 
the  results  compared  with  the  original  values.   No  measurable  super- 
saturation  effect  was  observed  in  these  experiments. 

In  order  to  investigate  the  relative  importance  of  surface 
effects  in  the  saturation  phenomena,  a  relation  between  the  solubility 
of  a  gaseous  solute  in  a  liquid  solvent  for  a  curved  surface  and  that 
for  a  planar  surface  was  derived.   The  complete  derivation  from  the 
appropriate  thermodynamic  relations  has  been  given  in  Appendix  1.   The 
ratio  of  the  solubility  on  the  curved  surface  to  that  on  the  planar 
surface  was  calculated  by  means  of  the  relation  obtained  for  the  0„  -  K  0 
system.   It  was  found  that  surface  effects  were  negligible  except  for 
bubble  radii  below  about  10   cm.   Consideration  of  the  numerical  value 
of  the  quantities  involved  in  the  relation  for  surface  effects 
suggests  that  the  error  will  be  approximately  the  same  or  smaller  for 
strong  KOH  solutions  and/or  higher  temperatures. 
Transfer  of  Samples: 

S:..-nples  of  the  saturated  solution  were  withdrawn  by  means  of  a 
gas-tight  Hamilton  syringe  ^.'ith  a  Luer  lock  to  which  was  attached  a 
Hamilton  teflon  valve  also  wich  a  Luer  lock.   While  withdrawing  samples 


DO 


from  the  saturating  vessel,  this  valve  fitted  on  the  valve  attached 
to  the  teflon  needle  dipping  in  the  solution.   The  samples  were 
withdrawn  extremely  slowly  so  that  at  no  stage  of  the  sampling  process 
was  the  flask  or  syringe  placed  under  reduced  pressure.   Several 
samples  v;ere  taken  and  rejected  before  the  final  sampling  was  made. 
As  a  check  on  the  procedure  of  sampling  and  degassing,  blank  samples 
consisting  of  KOH  solution  saturated  with  carrier  gas  were  injected. 
No  peak  was  obtained,  indicating  that  negligible  contamination  with 
atmospheric  or  other  gases  occurred. 


o 

For  experiments  at  temperatures  of  80  C  and  below,  two  pre- 

saturators  containing  the  KOH  solution  studied  were  used.   For  the 
experiments  at  100  C  additional  presaturators  were  necessary.   To  check 
that  the  solute  gas  stream  contacting  the  sample  was  fully  presaturated, 
samples  of  the  exit  gas  were  analyzed  for  water  vapor.   In  addition, 
the  concentration  of  KOK  in  the  solution  contained  in  the  saturation 
flask  was  checked  at  the  end  of  the  experiment. 

At  temperatures  of  100  C,  concentrated  KOH  solutions  rapidly 
attacked  glass  frits.   At  this  temperature  teflon  frits  and  teflon 
presaturators  were  used.   The  whole  assembly,  consisting  of  presatura- 
tors and  the  saturating  vessel,  was  completely  immersed  in  an  insulated 
stainless  steel  water  or  oil  bath.   The  temperature  of  the  bath  was 
controlled  to  -  0.05  C  by  means  of  a  proportional  temperature  controller, 
Thermometers  checked  for  accur: cy  against  NBS  calibrated  thermometers 
were  used  for  measuring  the  temperature. 
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Method  I  was  used  to  measure  the  solubility  for  systems  where 

this  v;as  above  about  10"'^  g.ir.ole  of  solute  gas  per  liter  of  eiectrolyce 

solution.   However,  the  solubility  values  for  SF  and  neopentane  in 

6 

electrolyte  solutions  having  more  than  207.,  by  weight  of  KOH  were  much 
less  than  10   ,  and  thus  could  not  be  measured  accurately  by  method  I. 
Therefore  a  modification  of  the  method  was  necessary.   This  consisted 
of  incorporating  a  concentration  technique  in  method  I  as  explained 
belovj. 
Method  II: 

The  principal  difficulty  in  extending  method  I  to  lesser 
solubilities  was  that  samples  larger  than  30-40  ml.  could  not  be 
injected  without  loss  of  accuracy,  because  of  broadening  and  tailing  of 
peaks.   This  difficulty  was  overcome  by  concentrating  the  stripped  gas 
by  means  of  lov;  temperature  adsorption  (58)  before  it  entered  the 
chromatograph.   This  provided  a  means  for  stripping  larger  volumes  of 
saturated  solution  of  the  dissolved  gas  and  permitted  the  rapid  intro- 
duction of  the  concentrated  gas  into  the  chromatograph. 

Figure  6  shows  a  schematic  diagram  of  the  apparatus  used  in 
method  II.   The  saturating  vessels  (not  shown),  the  stripping  cell  and 
the  drying  columns  were  merely  enlarged  versions  of  those  described 
in  method  I,   Solenoid  valves  1  and  2  were  set  to  perm.it  helium  carrier 
gas  CO  flow  directly  through  to  the  sample  column  of  the  chrom^atograph. 
Solenoid  valves  3  and  4  were  set  to  permit  stripping  helium  from  the 
drying  columns  to  flow  through  a  concentrator  and  a  moisture  trap  and 
then  to  the  atmosphere.   The  concentrator  was  a  1/4"  dia. ,  U-shaped 
copper  tubing  containing  a  6"  length  packed  with  activated  charcoal  in 
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its  central  portion  to  adsorb  the  solute  gas  at  liquid  nitrogen  tempera- 
ture, while  allowing  the  stripping  gas  to  pass  without  adsorption. 

Stripping  helium  was  first  purified  by  an  ETI  Standard  Laboratory 
Model  helium  diffusion  cell  and  then  passed  through  a  6'  long  1/4" 
dia.  copper  tubing  kept  in  a  liquid  nitrogen  bath  to  remove  all  adsorbable 
impurities.   It  was  then  passed  through  the  stripping  cell  and  the 
drying  column.   The  stripping  gas  together  v/ith  the  stripped  solute  gas 
flowed  through  solenoid  valve  3  to  the  concentrator  which  had  been 
immersed  in  liquid  nitrogen.   The  stripped  gas  was  adsorbed  in  the 
concentrator,  while  the  stripping  helium  passed  on  through  solenoid 
valve  4  and  the  moisture  trap  from  which  it  was  vented  to  the  atmosphere. 
The  moisture  trap  which  was  a  test-tube  immersed  in  a  liquid  nitrogen 
bath  was  necessary  to  guard  against  the  possibility  of  diffusion  of 
atmospheric  moisture  into  the  concentrator. 

The  saturated  KOH  solution  from  the  saturator  was  injected 
in  any  desired  a...ount  into  the  stripping  cell  and  the  solution  allowed 
to  strip  for  the  desired  time  (usually  about  5  min.),  the  concentrator 
having  been  previously  cooled  to  liquid  nitrogen  temperature.   At  the 
end  of  the  stripping  process,  valves  3  and  4  were  set  to  stop  the  flow 
of  the  stripping  helium  and  to  connect  to  valves  1  and  2.   The  liquid 
nitrogen  bath  was  removed  from  the  concentrator  and  replaced  by  a  bath 
of  boiling  water  for  a  specified  period  of  time  (5-10  min.),  which  caused 
essentially  complete  desorption  of  the  adsorbed  gases.   Valves  1  and  2 
v;ere  then  set  to  direct  the  flow  of  carrier  helium  through  the  concentra- 
tor and  the  dissolved  gas  sample  was  sv7ept  into  the  chromatograph  for 
analysis.   The  calibration  factor  was  obtained  as  in  method  I  by 


TO 


following  che  above  p-ocedure  and  using  i  standard  volume  of  water 
saturated  with  the  solute  gas  instead  of  XOH  solution. 

It  may  be  noted  uhat  the  length  of  the  packing  in  the  concen- 
trator and  the  chromatographic  column  size  should  be  carefully  chosen 
so  as  to  avoid  base  line  instability  which  may  otherwise  occur  because 
of  sudden  fluctuations  in  the  carrier  gas  flow  rate  on  switching  the  gas- 
stream  to  pass  through  the  concentrator.   The  packing  in  the  concentra- 
tor must  be  suitably  chosen  so  that  it  adsorbs  the  solute  gas  completely 
and  reversibly.   As  a  check  on  the  adsorption  process,  the  following 
experiment  was  made.   A  sample  of  pure  solute  gas  was  injected  into  the 
stripping  cell  while  the  helium  carrier  gas  was  allowed  to  sweep  it 
into  the  chromatograph  through  the  concentrator  without  being  adsorbed. 
The  same  quantity  of  the  gcs   sample  was  again  injected  while  the  con- 
centrator had  been  immersed  in  the  liquid  nitrogen  bath.   No  peak  was 
recorded.   Carrier  helium  to  the  concentrator  was  then  stopped  and 
directed  directly  to  the  chromatograph.   The  solute  gas  in  the  concen- 
trator was  desorbed  by  the  procedure  indicated  previously  and  swept  into 
the  chromatograph.   A  peak  was  obtained  whose  area  was  equal  to  that 
obtained  without  adsorption  within  the  limits  of  experimental  error. 
This  showed  that  adsorption  of  the  solute  gas  was  complete  and 
reversible. 

Since  the  method  is  based  on  low  temperature  adsorption  it 
cannot  be  used  for  gases  which  boil  at  temperatures  lower  than  that 
of  carrier  gas.   Extremely  pure  stripping  gas  is  necessary  in  these 
experiments  because  very  slight  concentrations  of  impurities  become 
concentrated  ir.  the  concentrator.   M.  hod  II  was  successfully  used  for 
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the  measuremenc  of  the  solubility  of  0  .   However,  an  error  in  the 

2 

saturation  procedures  was  discovered  after  these  measurements  were 

made.   Therefore,  they  vjere  repeated  with  method  I. 

The  above  procedure  was  not  found  to  be  satisfactory  for  S? 

5 

or  neopentane  because  the  activated  charcoal  used  as  a  packing  material 

in  the  concentrator  adsorbed  these  gases  irreversibly.   During  the 

course  of  a  search  for  new  packing  material,  it  was  found  that  SF  and 

6 

neopentane  could  be  trapped  completely  and  reversibly  at  liquid 
nitrogen  temperature  with  a  colunm  containing  no  packing.   As  a  result 
no  heating  of  the  concentrator  was  necessary  to  recover  the  trapped  solute 
gas.   Thus  method  III,  described  below,  could  be  used  for  these  measure- 
ments . 
Method  III: 

This  method  is  very  similar  to  method  I.   A  schematic  diagram 
is  shown  in  Figure  7.   The  stripped  gas,  after  passing  through  the  drying 
columns,  was  directed  to  the  concentrator,  which  was  very  similar  to  that 
in  method  II  but  did  not  contain  any  packing.   The  gas  stream  from  the 
concentrator  passed  through  the  stainless-steel  coil  immersed  in  a 
water  bath  and  returned  to  the  chromatograph.   The  following  procedure 
was  used  to  make  a  measurement. 

KOK  solution,  saturated  with  the  solute  gas,  was  injected  into 
the  stripping  cell,  the  concentrator  having  been  previously  cooled  to 
liquid  nitrogen  temperature.   At  the  end  of  the  stripping  process 
(usually  5-10  min.)  the  liquid  nitrogen  bath  was  removed,  and  the 
condensed  solute  gas  rapidly  evaporated  and  passed  through  the  chroma- 
tograph to  give  a  sharp  peak. 
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There  was  a  slight  base  line  disturbance  when  cooling  of  the 
concentrator  by  means  of  liquid  nitrogen  was  started  and  again  when 
the  liquid  nitrogen  bath  was  removed;  however  these  disturbances  died 
cut  almost  instantaneously  and  did  not  affect  the  results. 
Calculation  Procedure: 

The  following  equation  was  derived  in  Chapter  2. 

P2  =  VZ^'Z  ^^^^^ 

Let     p'  =  vapor  pressure  of  water  at  25  C 

p  =  vapor  pressure  of  KOH  solution  at  T°C 

o         o 
X  (P)  and  x^(P-p')  =  mole  fraction  of  solute  gas  in  water  at  P  and  P-p 

solute  gas  pressure,  respectively. 

X  (?)  and  X  (P-p)  =  mole  fraction  of  solute  gas  in  KOH  solution  at  P  and 

P-p  solute  gas  pressure,  respectively. 


o  o 


P  =  K272^2^''^  ^^^^^ 


or 


x°(P)  =  -^  (155) 

V2 

For  fairly  low  temperatures  where  p  and  p'  are  small  K  can  be 

considered  constant  with  respect  to  pressure  and 


x°(P-p')  =^  (156) 

K272 

or 


4^-p-p')      p.v 


x°(P) 


(157) 


By  a  similar  assumption  regarding  K  in  case  of  KOH  solution, 


From  (157)  and  (Ipcl) 

x°(P-p')      ,   x^CP)   A° 
x^CP-p)   P-p    X2(?)   A2  ^^^ 

Where  A„  and  A  are  peak  areas  for  water  and  KOH  solution  respectively, 
From  (159) 

X2(P)=x°(P)(^^  Iff  (160) 


CHAPTER  k 
RESULTS 

The  experimental  data  for  the  solubility  of  various  gases  are 
presented  in  Table  1,  where  activity  coefficients  of  the  solute  gases 
in  KOH  solutions  have  been  tabulated.   These  values  have  been  calculated 
from  the  experimental  results  of  the  solubility  of  a  given  gas  in  water 
and  in  KOH  solutions  by  means  of  the  following  equation  obtained  from 
equations  (22)  and  (23)  of  Chapter  2. 

x^(?)  is  the  mole  frc.ction  of  the  solute  gas  dissolved  in  a 
KOH  solution  in  equilibrium  with  the  solute  gas  at  I  atmosphere  partial 
pressure,  and  is  calculated  from  the  experimental  measurements  by 

using  the  following  equation  derived  in  Chapter  3- 

>'2W-K7)(K^) 

Vapor  pressure  data  for  KOH  solutions  were  taken  from  International 
Critical  Tables  (59) . 

In  all  cases,  the  solubility  values  are  the  mean  of  four  or 
more  replicate  measurements.   Tha  absolute  accuracy  of  the  reported 
data  depends  upon  the  amount  of  the  dissolved  gas  present  in  a  fixed 
quantity  of  electrolyte  solution  which,  for  a  given  gas,  depends  upon 
the  electrolyte  concentration.   Considering  the  errors  involved  in 
saturation  of  the  electrolyte  samples,  transferring  of  the  samples, 
chromatographic  analysis  and  KOH  concentration  measurements,  the 
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overall  accuracy  of  the  reported  results  is  estimated  to  be  about  -  25© 

for  solute  gas  concentrations  up  to  10   g. mole/liter ,  and  about  -  6^ 

q      _  ^ 
in  the  concentration  range  of  10   to  10   g. mole/liter .   (See  Appendix 

-3 
2.)   Solubility  values  for  pure  water  lie  in  the  range  or  10   to 

-L 
10  '  g. mole/liter  and  the  results  reported  here  in  most  cases  agree 

with  the  literature  values  within  -  2^p.      This  supports  the  accuracy 

quoted  above  for  the  solubilities  in  dilute  electrolyte  solutions. 

The  estimate  of  the  accuracy  for  the  concentrated  electrolyte  solutions 

(dissolved  gas  'v  10   -  10   g. mole/liter)  is  difficult  to  confirm 

because  few  reliable  data  are  available.   However,  the  agreement  between 

the  solubilities  of  oxygen  at  25  C  obtained  in  this  work  and  those  of 

Davis  et_  al.  (60)  and  Geffcksn  (6l)  (discussed  later)  suggest  that  the 

accuracy  estimations  are  reasonable. 

The  range  of  temperature  and  potassium  hydroxide  concentration 
over  which  the  solubility  of  each  gas  was  studied  is  shown  in  Table  2. 

The  activity  coefficients  of  che  solute  gas  at  each  temperature 
have  been  plotted  on  a  semilog  paper  w:_th  abscissa  as  potassium 
hydroxide  concentration  in  g. mole/ liter ,  and  Figures  8-l4  show  these 
plots.   The  semilog  plot  has  been  used  since  most  of  the  electrostatic 
theories  of  saiting-out  (Chapter  2)  predict  a  linear  behavior  on  such 
a  plot  at  moderate  electrolyte  concentrations. 
Comparison  of  Results: 

Measurements  of  the  solubility  of  oxygen  and  hydrogen  in  aqueous 
potassium  hydroxide  solutions  have  been  made  by  Geffcken  (6l)  at  15 
and  25  0  for  the  concentracion  range  of  0  to  l.h   molar,  and  by  Knaster 
and  Apel*baum  (62)  ac  21,  ^5  and  75  0  for  concentrations  up  to  10  molar. 
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TABLE  2 
EXPERIMENTAL  SYSTEMS  AND  CONDITIONS 


/^ 

KOK 

Concentration 

Solute  Gas 

Tianperature 

,  °c 

Range,  1, 

/t .  ^' 

...  ■  ..i  .       -  — 

Figure  No. 

°2 

25, 
100 

40, 

DO, 

8c 

0  - 

ko 

■  50 
-  56 

8 

8 

«2 

25, 
100 

40, 

6o, 

80 

0 

38 

-  50 

-  56 

9 
9 

He 

25. 

40. 

60, 

80 

0 

-  30 

10 

Ar 

25. 

40, 

60, 

80 

0 

-  40 

11 

S^6 

25. 

40. 

6c, 

80 

0 

-  30 

12 

CHj, 

25, 

40. 

60, 

80 

0 

-  40 

13 

neo- 

•C.H^o 

25. 

40. 

60, 

80 

0 

-  30 

14 
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Figure  8.   Activity  Coefficients  of  Oxygen. 
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Figure  9-   Activity  Coefficients  of  Hydrogen 
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Figure  11.   Activity  Coefficients  of  Argon. 
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Figure  12.   Activity  Coefficients  of  Sulfur 
Hexaf luoride . 


36 


iog'_ 


^^"^ 


•r-l 

o 


a 
o 
u 


o 


IC 


o 

25   C 

A 

itO 

a 

60 

o 

80 

2^+68 
KOH  Concentration,  G. Mole/Liter 
Figure  13-   Activity  Coefficients  of  Methane, 


10 


8? 


lOG 


?-^ 


c 
aj 

•r-l 

o 


0 
o 


•r-l 
> 

O 


10  __ 


- 

0    25°C 
A    ^0 

__ 

Q    '° 

0    80 

1 

0        2        4        6        8 

KOH  Concentration,  G. Mole/Liter 
Figure  lii-.   Activity  Coefficients  of  Neooentane. 


88 


Davis  e_C  a  I .  (60)  have  reported  oxygen  solubility  measurements  at 
temperatures  of  0,  25  and  6C  C  for  the  concentration  range  of  G  to  12 
g. mole/liter.   More  recently,  Ruetschi  (21)  measured  the  solubilities 
of  hydrogen  in  KOH  solutions  at  30  C. 

To  facilitate  the  comparison  of  the  results  of  the  present 
investigation  with  those  of  the  previous  workers,  the  data  for  oxygen 
and  hydrogen  at  25  C  have  been  plotted  in  Figures  15  and  l6.   Also 
the  results  obtained  for  the  solubility  of  different  gases  in  water, 
together  with  the  literature  values  have  been  tabulated  in  Table  3. 
The  ^  difference  between  the  values  obtained  in  this  work  and  the 
literature  values  has  been  shown.   It  is  seen  from  this  table  that  the 
agreement  between  the  solubilities  reported  here  and  the  literature 
values  is  in  most  cases  better  than  2^. 
Solubilicy  of  Oxygen: 

Comparison  of  the  results  for  oxygen  at  25  C  with  those  of 
other  investigators  (Figure  15)  reveals  that  the  solubility  values 
obtained  in  this  work  agree  within  2^  and  p/o  with  those  obtained  by 
Geffcken  (6l)  and  Davis  et_   al.  (6C),  respectively.   The  experimental 
results  of  these  investigators  are  distributed  on  both  sides  of  the 
solid  line  of  Figure  15 .   The  agreement  with  the  results  of  Knaster 
and  Apel'baum  (62)  is  quite  good  at  very  low  and  very  high  KOH  con- 
centrations.  For  example,  for  KOH  concentration  of  1  molar  the  agree- 
ment is  better  than  3^.  whereas  for  10  molar  KOH  solution  it  is  about 
IC^.   Hovjever,  at  intermediate  concentrations,  the  disagreement  is  quite 
marked,  the  results  of  Knaster  and  Apel'baum  falling  consistently 
belov;  the  solid  line.   At  temperatures  other  than  25  C,  the  values 
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reported  by  various  workers  differ  appreciably.   The  disagreement  is 
quite  pronounced  even  at  low  KOH  concentrations.   For  example,  the 
solubility  reported  by  Davis  ez_   3I .  (cO)  in  pure  water  at  cO  C  differs 
significantly  (about  2C^)  from  the  currently  accepted  literature  values 
(63),  and  this  suggests  that  discrepancies  r::ay  be  present  at  other 

concentrations  also.   The  agreement  wich  the  results  of  Knaster  and 

>.  o  o 

Apel'baum  at  ^rO  and  6C  C  is  of  the  saae  general  nature  as  at  25  C. 

From  Table  3»  ^-^    is  seen  that  the  agreement  for  the  results  in  pure 
water  is  within  2'p. 
Solubility  of  Hydrogen: 

Comparison  of  the  results  for  hydrogen  at  25  C  shows  (Figure 
16)  that  solubilities  obtained  in  this  work  are  in  fairly  good  agree- 
ment with  those  of  other  workers.   The  agreement  with  the  results  of 
Geffcken  (6I)  is  within  2^  whereas  with  those  of  Knaster  and  Apel'baum 
(62)  is  within  lOp.   For  dilute  KOH  solutions,  the  reported  results 
agree  quite  well  (within  S'/o)  with  those  of  Ruetschi  (21).   However,  the 
values  of  activity  coefficients  obtained  by  him  are  consistently 
lower  than  those  anticipated  by  previous  workers  and  the  deviations 
between  his  values  and  those  obtained  in  this  work  at  higher  KOH 
concentrations  are  of  the  order  of  IC73.   At  temperatures  of  k-O   and 
6G  C,  the  disagreement  between  the  reported  results  and  those  of  Knaster 
and  Apel'baum  (6I)  is  greater.   The  agreement  between  the  solubilities 
for  pure  water  obtained  in  this  work  and  the  literature  values  is  in 
most  cases  vjithin  2^. 
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Figure  16 .   Activity  Coefficients  of  Hydrogen  at  25  C. 
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Solubility  of  Helium  and  Argon: 

No  data  for  the  solubility  of  argon  and  helium  in  aqueous 
potassium  hydroxide  solutions  are  available.   However,  a  considerable 
amount  of  work  has  been  done  on  the  solubilities  of  these  gases  in 
water.   Table  3  shows  that  results  reported  here  for  the  solubility 
of  helium  in  V7ater  are  in  excellent  agreement  (better  than  1.6^) 
with  those  obtained  by  Morrison  and  Johnstone  (64)  and  Lannung  (65) • 
The  agreement  between  reported  values  of  solubility  for  argon  and 
those  of  Morrison  and  Johnstone  (o^r)  is  of  the  order  of  2.5^- 
Solubility  of  Sulphur  Kexaf louride : 

As  in  the  case  of  helium  and  argon,  no  data  are  available  for 
the  solubility  of  this  gas  in  potassium  hydroxide  solutions. 
Friedman  (06)  has  reported  the  solubility  of  this  gas  in  water  at 
2o   C  (see  Table  3)  which  agrees  with  the  corresponding  value  reported 
here  within  Vjo. 
Solubility  of  Methane: 

As  with  the  last  three  cases  comparison  can  only  be  made  with 
the  solubility  of  this  gas  in  water.   The  results  obtained  here  at 
various  temperatures  are  in  excellent  agreement  (wichin  1.4^)  with 
those  commonly  reported  in  the  literature  (63) • 
Solubility  of  Neopentane: 

Wetlaufer  et_  al .  (67)  have  reported  the  solubility  of  neopen- 
tane in  water  at  Ip ,  35  and  i+5  C.   The  agreement  between  these  values 
and  those  obtained  in  this  work  is  poor.   However,  the  solubility  values 
reported  by  these  authors  when  plotted  against  temperature  do  not 
yield  a  smooth  curve. 


CHAPTER  5 
DISCUSSION 

The  results  of  Table  1  reveal  that  the  activity  coefficients 
in  aqueous  potassium  hydroxide  solutions  of  all  the  gases  studied  in 
this  v7ork  increase  very  rapidly  with  increase  in  potassium  hydroxide 
concentration.   Owing  to  the  reciprocal  relationship  between  the 
activity  coefficient  and  the  solubility,  the  above  statement  implies 
that  solubility  decreases  very  rapidly  with  increase  in  potassium 
hydroxide  concentration.   As  an  illustration,  the  solubility  of 
oxygen  in  pC/S  by  weight  aqueous  potassium  hydroxide  solution  is  less 
than  I'p   of  the  value  in  v;ater.   The  salting-out  effect  is  greatest  for 
sulfur  hexafluoride  and  neopentane  and  least  for  helium  and  hydrogen. 
A  comparison  of  the  salting-out  of  some  of  these  gases  in  various 
electrolytes  (68)  indicates  that  this  effect  is  more  pronounced  for 
potassium  hydroxide  than  for  other  electrolytes.   From  Figures  Q-lk   it 
may  be  seen  that  the  solubilities  of  these  gases  change  with  tempera- 
ture.  The  temperature  effect  is  most  marked  for  dilute  potassium 
hydroxide  solutions,  when  the  solubility  decreases  with  increasing 
temperature,  whereas  for  concentrated  solutions  the  solubility  either 
remains  constant  or  increases  slightly. 
5  . 1 Classical  Thermodynamics 

In  section  2.1,  phase-equilibrium  relations  for  a  gas  dissolved 
in  an  aqueous  electrolytic  solution  were  derived.   These  relations  are 
useful  in  relating  the  solubility  phenomena  to  the  physical  characteris- 
tics of  the  solute  gas  and  the  concentration  of  the  electrolyte. 

9k 
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Partial  Molal  Heats  of  Solutioa 

The  effect  of  temperature  on  solubility  depends  upon  the  con- 
centration of  the  electrolyte  and  the  type  of  solute  gas.   This  effect 
can  be  quantitatively  examined  by  calculation  of  partial  molal  heats 
of  solution. 

The  following  relation  was  derived  for  the  partial  molal  heat 
of  solution  in  Chapter  2. 

,  ain(7pK  )  .         m 
^     °^    ^P,n.     RT^ 

The  symbols  have  been  defined  previously.   The  quantity 
Z!^I  =  H  -  h^  is  the  change  in  enthalpy  when  1  mole  of  gas  in  the 
ideal  gaseous  state  dissolves.   AH^  values  for  various  gases  were 
obtained  by  plotting  the  product  (7  K  )  against  reciprocal  temperature 
on  semilog  co-ordinates,  and  using  equation  (32).   The  calculated 
values  of  Z\H,p  are  tabulated  in  Table  h.      A   typical  plot  of  partial 
molal  heats  of  solution  is  that  for  oxygen  and  has  been  shown  in 
Figure  I7.   It  is  seen  that  the  heat  of  solution  is  negative  at  lower 
potassium  hydroxide  concentrations  and  becomes  less  negative  as  the 
pocassium  hydroxide  concentration  increases.   Also,  at  lower  tempera- 
tures the  heat  of  solution  is  more  negative.   Similar  behavior  vjas 
exhibited  by  other  gases  studied  in  this  work  and  such  plots  for  these 
gases  are  included  in  Appendix  3  (Figures  27-32) . 

As  seen  from  Table  4,  for  a  given  temperature  and  potassium 
hydroxide  concentration,  Zj>-„  varies  with  the  solute  gas.   It  is  reasonable 
to  expect  that  this  variation  of  partial  molal  heats  is  due  to  physical 
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characteristics  of  the  solute  gas,  or  at  least,  to  its  effect  on  the 
properties  of  the  solvent.   To  this  end  zXH,^  was  plotted  against  e/k  in 
Figure  l8,  £  being  the  energy  parameter  (kd)    of  the  Lennard- Jones 
6,12  potential.   /VH  seems  to  correlate  well  with  this  parameter.   It 
is  seen  that  Ui^   decreases  as  e/k  increases  for  lower  concentrations 
of  KOH  but  exhibits  a  minimum  at  higher  concentrations.   It  may  be 
noted  that  ZJi^  was  also  plotted  against  the  Lennard-Jones  C   and  the 
polarizability,  p  ,  and  that  satisfactory  correlations  were  also 
obtained  using  these  force  constants. 
Partial  Molal  Energies  of  Solution 

It  has  been  shown  (39>  ^)  that  the  energy  of  solution  of  a 
gas  molecule  in  a  solvent  is  an  important  variable  in  correlating  and 
predicting  gas  solubilities.   This  quantity  is  given  by  equation  (kO) 
of  Chapter  2  which  is 

The  above  relation  was  obtained  by  neglecting  the  partial  molal  volume 

-L  g  , 

or  the  dissolved  gas  V^  as  compared  to  the  volume  of  the  gas,   V^  (since 

— L  <^<;;  V,,)  ,  and  assuming  that  the  gas  phase  behaves  ideally  at  the  low 

2     '^ 
pressures  encountered  in  this  work.   The  calculated  values  of  AE  have 

been  tabulated  in  Table  5-   Energies  of  solution  for  all  the  gases  have 

been  plotted  against  temperature;  Figure  I9  shows  such  a  plot  for  0    . 

It  is  seen  that  Z\E  is  negative  at  low  potassiiam  hydroxide  concentrations 

and  temperatures,  and  that  it  tends  to  increase  as  the  potassium  hydroxide 

concentration  or  temperature  increases.   The  plots  of  energies  of 

solution  for  other  gases  show  a  similar  behavior,  and  are  presented  in 
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Appendix  3  (Figures  33-38)-   To  examine  the  variation  of  energy  of 
solution  with  the  solute  gas,  this  quantity  was  plotted  against  e/k 
in  Figure  20.   It  appears  that  there  is  a  definite  correlation  between 
AEp  and  e/k.   As  in  the  case  of  /^H  ,  /^  decreases  with  increase  in 
e/k  for  the  solubility  of  these  gases  in  water,  but  it  exhibits  a 
minimum  at  high  concentrations. 
Partial  Molal  Entropies  of  Solution 

The  following  relation  between  AS  ,  the  partial  molal  entropy 
of  solution  and  Z^  ,  the  partial  molal  enthalpy  of  solution  was 


obtained  in  Chapter  2. 


-    ^2 


2    T 

The  calculated  values  of  AS-  are  tabulated  in  Table  6.   A  typical  plot 
of  ASp  as  a  function  of  temperature  is  shown  in  Figure  21  for  0  .   It 
may  be  seen  that  the  entropy  change  is  negative  at  lower  potassium 
hydroxide  concentrations  and  becomes  less  negative  as  potassium 
hydroxide  concentration  increases.   This  implies  that  the  system  becomes 
more  organized  after  dissolution  of  the  gas  c.nd  that  it  becomes  less 
organized  as  the  concentration  of  potassium  hydroxide  increases . 
Dependence  of  AS^  on  temperature  and  KOH  concentration  for  other  gases 
was  very  similar  to  that  of  oxygen  and  is  shown  in  Appendix  3  (Figures 
39-^4-^1).   To  show  the  effect  of  the  solute  gas  on  ASp,  this  quantity 
was  plotted  against  e/k  (Figure  22).   Similar  variations  were  found  as 
for  the  quantities  AKp  and  Z\E  . 

5 .2 Test  of  Electrostatic  Theories 

In  Chapter  2,  a  summary  was  given  of  the  electrostatic  theories 
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which  have  been  proposed  by  various  investigators.   In  this  section, 
some  of  these  theories  are  used  to  perform  calculations  for  the  systems 
studied  in  this  work. 
5.2.1   Theory  of  Dabye-McAulay 

S°  -   s   n  e^ 

2  o  1     J 

D-D 

All  symbols  have  been  defined  previously,   p  =  — is  an 

o 
empirical  constant  which  should  be  obtained  experimentally.   However, 

these  data  are  not  often  available  and,  since  this  situation  arises 

in  alraost  all  the  electrostatic  theories,  we  resort  to  theoretical 

expressions  which  predict  the  effect  of  concentrations  of  the  solute 

molecule  on  the  dielectric  constant  of  the  mixture  of  the  solute  and 

water. 

In  the  case  of  highly  polar  systems,  Kirkwood's  (20)  limiting 

equation  of  dielectric  constant  is 


K 

where   n  =  number  of  molecules  of  component  k/volume. 

jj.^   =  moment  of  the  individual  noneleccrolyte  molecule 

^  =  total  moment  of  the  central  molecule  and  its  neighbors 

in  a  spherical  region  within  which  the  dielectric  constant 
differs  significancly  from  its  macroscopic  value. 

a  ,  =  polarisation  resulting  from  a  displacemenc  of  charge 
within  the  solute  molecule. 

By  following  the  treatment  of  Edsalland  Wyman  (20)  it  can  be  shown 

that 


108 


D   -  D   V      P 

J^^^^^'l^  (163) 


o 

2     =   total  molar  polarisation  of  the  nonelectrolyte  and  is 
defined  by  the  following  equation: 

where  a     is  the  polarisability  of  the  solute  gas  molecule. 

Substituting  P   in  equation  (l63)>  the  following  relation  for  p  for  a 

nonpolar  gas  is  obtained  (^  =  O) : 

D-D   V^   6tkx 
J  =  ^ =  —  --Tr^  (165) 

^    D  n     '^     D  \   ^z 

o      N      o 

In  equation  (^9)  n.  is  number  of  ions  of  type  j  per  cm.   of 
solution  and  is  given  by: 

V.C.N 

n.  =  -tL-J-  (166) 

J    1000  ^    ^ 

where   C  -  concentration  of  electrolyte  in  g. moles  per  litre. 

V.    =   number  of  ions  of  type  j  per  molecule  of  electrolyte. 

Thus  the  equation  of  Debye  and  McAulay  can  be  written  as 

2 
8  s     V .  e  . 

^°Sig'^2  =   ^   2.303x2CC0kTD        f    "b"^  ^   S  ^^^^^ 

o        J 

=  Vs 

2 

-        V  .e  . 

"^^^^  '^S   =  2.303x20C0kTD^  -  "V  ^^^^^ 

The  above  equation  was  used  to  calculate  the  constant  k  ,  the 

o 

salting-out  coefficient,  for  all  the  gases  studied  in  this  work  at  25  C. 

The  calculated  results  are  tabulated  in  Table  7.   Values  of  b.  used  are 

J 
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Gas 
He 

\ 
Ar 

°2 

CH, 


^^6 
neo-C_H 


TABLE  7 


SALTING-OUT 

COEFFICIENT, 

kc  at  2p  C 

Debye- 

-McAulav 

Conway 

A 

B 

ExDt . 

0.0173 

0.035 

0.081 

0.0?? 

O.CirS 

O.C89 

0.09l{- 

0.129 

0.05^^ 

0.111 

0.099 

0.179 

0.072 

O.II.7 

0.108 

0.180 

0.085 

0.17^+ 

0.112 

0.197 

0.23.5 

0.i;32 

0.187 

0.33 

- 

- 

- 

0.3i^ 

5  12 


A  -  Results  calculated  using  crystal  radius  of  the  ion 

B  -  Results  calculated  using  ionic  radius  calculated  from  ionic 
mobility  data 
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given  in  Table  8.   Since,  the  theories  of  Bulter  and  Belton  give 
expressions  v.'hich  are  very  similar  to  that  of  Debye-McAulay ,  they  v;ere 
not  used  to  predict  salting-out  coefficients. 
3.2.2       Theory  of  Conway  et  al. 

Conway's  theory  (see  Chapter  2)  gives  the  following  expression 
for  the  salting-out  coefficient  for  a  nonpolar  gas  in  1:1  electrolyte: 

S°  -  S    no^n-,  +  n^^   Z?£^ 


S°  C 
^2  ^S 


^  18("1  "•  "2)    2     ^2        r-  /  1 1 N 

ICOCd  -  C^W  ^   200GkTD   _  \   (l)  "*"   (2)  / 
S  o      ^h     ^h 

-(I)"     (169) 

where  n.  and  r^   are  the  primary  hydration  number  and  radius 
respectively  for  the  ion  i,  R  was  defined  in  Chapter  2. 

S°  -  S^    Iny, 

Lim     -r ^  =  — pr"  =  2.303k„  (170) 

C^-.C  S°C3 


or 


1    r  'Q("i  ^  '-^  .  ^2^  ^2     /I      .1      M 


,2  2r, 

L(i)  ^  ( 


S   2.303  I     lOOCd       2000kTD   V  (l)    (2)  7  . 

°   ^h     "h     " 

(171) 

+      — 
The  hydration  radii  and  hydration  numbers  for  k  and  OH  ions 

used  to  calculate  the  salting  coefficients  for  various  solute  gases  are 

listed  in  Table  8.   The  hydration  number  for  OH  ions  was  not  available 

in  the  literature  and  was  calculated  from  the  partial  molal  entropy  of 

hydration.   Ulich  (7I)  (see  also  reference  72)  assumed  that  the  partial 

molal  entropy  of  hydration  arises  froni  the  restricted  freedom  of 
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TABLE  8 


VALUES  OF  PARAMETERS  b,  r,  AND  n 

n 


Ion 

b!^)A° 

b(^)A= 

,o 

r^_,A 

n 

K^ 

1.33 

1.25 

2.72^"^ 

k.iM 

oh" 

1.76 

OAc 

3.00^^^ 

5.3<^' 

(a)  Crystal  radii  (references  69  and  7C) 

(b)  Radii  calculated  from  ionic  mobility  (reference  70) 

(c)  Reference  19,  (d)   Reference  70  (e)   Calculated  in  this  work 

TABLE  9 
PARTIAL  MOLAL  VOLUMES  OF  GASES  IN  WATER  AT  25°C 

Gas 

He 

«2 
Ar 

°2 

neo-C^H^2 

* 

Estimated 


V2,cm.3/g 

.mole 

Ref 

erence 

7.7 

19 

19.4 

73 

24.2 

19 

32.0 

Ih 

36.0 

75 

105* 

112 


movement  of  only  those  v;ater  molecules  contained  in  the  primary  shell 
of  hydration.   The  lowering  of  the  entropy  per  mole  of  water  is  assumed 
equal  to  the  loss  of  entropy  accompanying  the  transformation  of  the 
liquid  to  a  solid,  or  to  a  decrease  of  entropy  during  the  formation  of 
solid  hydrates  (-6  to  -7  cal/de-) ,  the  primary  hydration  numbers  are 

given  by :  _ 

^^h 
n=  --^ 

where   AS^  =  S   -  ? 
h    aq    g 

S„  =  partial  molal  entropy  of  the  ion  in  solution 

S   =  partial  molal  entropy  of  the  ion  in  gaseous  state 

The  term  S^,  for  an  ion  which  is  monatomic  in  the  vapor  state  at  296°K 

o 

can  be  calculated  from  the  Sacker-Tetrode  equation  (72): 

S^  =  26.0  +  6. So  LogM  (172) 

M  =  Molecular  weight 

The  partial  molal  entropy  of  the  OH  ion  in  solution,  S   ,  was  obtained 

aq 

from  the  literature  (69).   Table  9  lists  the  partial  molal  volume  of 
different  gases  in  water  used  to  calculate  the  results  by  equations 
(168)  and  (171). 

The  salting-out  coefficients  calculated  by  Conway's  theory  are 
presented  in  Table  7-   For  comparison  this  table  also  shows  the 
experimental  results  and  the  values  obtained  by  the  Debye-McAulay 
theory.   The  results  for  neopentane  could  not  be  calculated  due  to 
lack  of  partial  molal  volume  data  and  therefore  are  not  included  in 
Table  7-   It  may  be  noted  that  the  theory  of  Debye-McAulay  predicts  a 
linear  relationship  between  the  logarithm  of  the  activity  coefficient 
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and  the  molar  concentracion  of  the  electrolyte.   Similar  behavior  is 
also  predicted  by  Conway's  theory  at  very  low  electrolyte  concentrations. 
The  plots  of  Figures  S-14  reveal  that  these  theories  agree  with  the 
experiment  in  so  far  as  linear  behavior  is  concerned;  however,  there  is 
a  considerable  disagreement  between  theoretical  and  experimental  values 
of  salting-out  coefficients.   The  theories  of  both  Debye-McAulay  and 
Conway  predict  salting-out  coefficients  which  are  much  lower  than  the 
experimental  values.   In  addition  the  variations  of  this  coefficient 
with  the  size  of  the  solute  molecule  are  not  correctly  predicted.   In 
the  case  of  the  Debye-McAulay  theory  such  discrepancies  were  recognized 
by  Randall  and  Failey  (10)  for  several  other  solute-electrolyte 
systems , 

In  the  Debye-McAulay  theory  b.  (equation  49)  is  the  radius  of 
the  ion  j.   However,  the  method  of  obtaining  b.  is  a  controversial  ques- 
tion (17,  19).   For  this  reason,  results  in  this  work  were  calculated 
by  using  both  crystal  radii  and  radii  calculated  from  ionic  mobility. 
In  Table  7,  it  is  seen  that  the  radii  calculated  from  ionic  mobility 
generally  give  better  results  than  the  crystal  radii.   In  fact  the 
Debye-McAulay  theory  in  most  cases  gives  better  results  than  Conway's, 
provided  radii  calculated  from  ionic  mobility  are  used.   The 
choice  of  the  hydration  number,  n  and  the  hydration  radius,  r   in 
the  case  of  Conway's  theory  is  as  difficult  as  was  the  case  for  b.  in 
the  Debye-McAulay  theory.   Since  these  parameters  are  difficult  to 
obtain  there  are  often  large  variations  in  the  values  reported  by 
various  investigators  (19,  70).   The  values  c  "     d  r,_  used  in  this 
work  are  those  which  were  believed  to  -^  . ^^reeraent  with  the 
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experimental  data. 

The  failure  of  electrostatic  theories  to  predict  the  salt 
effects  correctly  can  be  attributed  to  the  simplifying  assiiraptions 
made  in  these  theories.   One  of  the  most  drastic  oversimplifications 
is  the  assumption  that  the  solvent  is  a  continuous  medium,  having 
constant  dielectric  properties.   In  reality,  the  solvent  is  a  system 
of  molecules,  and  the  interaction  energy  of  the  system  must  be  obtained 
by  considering  individual  molecular  interactions.   The  electrostatic 
theories  of  Debye-McAulay  (7)>  Debye  (ll)  and  Conway  (I9)  attempt  to 
explain  the  salt  effects  in  terms  of  electrical  energy  stored  in  an 
electrolyte  system  and,  in  fact,  appear  to  lay  too  much  emphasis  on 
this  energy.   However,  the  quantity  of  interest  is  the  change  in  this 
energy  due  to  the  addition  of  the  solute  molecule,  and  this  change  is 
assumed  to  be  given  solely  by  the  change  in  dielectric  constant  of  the 
solvent  which  takes  place  when  the  solute  molecule  is  dissolved.   In 
other  words,  the  change  in  the  free  energy  of  the  system  which  takes 
place  due  to  the  molecular  interactions  of  the  solute  and  the  solvent 
species  is  never  implicity  calculated.   Consequently,  since  the  gases 
have  a  generally  small  effect  on  the  dielectric  constant  of  the  solvent 
(21),  but  interact  strongly  with  the  solvent  species,  they  are 
predicted  to  show  erroneous  salting-out  coefficients. 

In  addition  to  the  above  objections  to  the  basic  model,  Debye 
(11)  and  Conway  ( I9)  assume  a  Boltzmann  distribution  for  the  presence 
of  solute  molecules  around  an  ion  which  at  best  is  only  an  approximation. 
The  exact  distribution  function  can  only  be  obtained  by  considering 
the  detailed  molecular  interactions.   In  such  circumstances,  further 
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attempts  to  improve  the  basic  theory  by  considering  dielectric  satura- 
tion effects  (by  introducing  the  concepts  of  hydration  shells)  as 
done  by  Conv/ay  ( I9)  seem  to  only  complicate  the  situation  and  can 
hardly  be  expected  to  yield  satisfactory  results. 
^  .3 Test  of  Scaled  Particle  Theory 

The  scaled  particle  theory  developed  in  section  2.2.ii-  is  used 
here  to  calculate  theoretical  values  for  the  solubility  of  various 
gases  in  potassium  hydroxide  solutions. 

The  following  equations  were  derived  in  Chapter  2 


^"(Va^=  fe^l*^"(f  S-j)  (128) 


_h  . 


Ll  -  C3    (1  .  C  )2^  ^  2    '  3  ^  kT  ^^   2  ^ 

(136) 


a.   1 


where   ^  =  ^  2  (  -^  )  p 
J 

p 

_h 
Evaluation  of  the  free  energy  of  making  the  cavity,  g   requires 

a  knowledge  of  hard  sphere  diameters  and  censity  of  the  solvent 

(potassium  hydroxide  solutions).   The  densities  of  potassium  hydroxide 

solutions  were  obtained  from  references  (76)  and  (77)-   The  hard  sphere 
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diameters  of  the  various  species  were  computed  as  shown  below. 
Determination  of  z   for  Water: 

The  value  of  a   for  v;ater  was  obtained  by  Pierotti  (37)  by 

_s 

the  follov7ing  procedure.   Since  g  is  related  to  the  poiarizability 

of  the  solute  molecules  (see  equation  (15£);  it  will  be  shown  later 

that  the  first  term  of  this  equation  is  also  related  to  polarizability) , 

a  plot  of  experimental  values  of  ln('y  K  )  against  the  polarizability, 

p  of  spherical  solutes  (the  rare  gases)  should  be  a  smooth  curve. 

Figure  23  shows  such  a  plot  for  the  solubility  of  rare  gases  in  water. 

It  was  shown  by  Pierotti  (36)  that  the  extrapolation  of  such  a  curve 

to  zero  polarizability  was  equivalent  to  obtaining  the  solubility  of 

a  hard  sphere  of  diameter  2.58  A  .   This  was  expressed  as: 

_h 

ln(^^K^)°  =  ||/ln(^Zn.)  (I73) 

p^-.0 

a   -*2.p8  A° 
a 

where  (-y  K  )   is  the  Henry's  law  constant  for  a  hypothetical  solute 
having  a  hard  sphere  diameter  of  2.58  A  .   Thus  O   for  water  was  cal- 
culated by  means  of  equation  (173)  using  a  value  of  2.58  A  for  O 
and  the  experimental  value  of  (-v  K  )   at  25  C.   This  value  is  shown 
in  Table  IC.   It  may  be  pointed  out  that  O   calculated  by  the  above 
procedure  will  vary  somewhat  with  temperature.   However,  initially 
results  will  be  calculated  using  the  constant  G  calculated  above  (2?  C) , 
and  the  temperature  dependance  of  a  will  be  discussed  later. 
c  for  Ions 

This  parameter  is  not  available  in  the  literature.   It  seems 
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TABLE  10 
PROPERTIES  AND  PARAMETERS  FOR  SOLVENT  SPECIES 


Solvent 

Parameter 

Polarizability 

Parameter 

Dipole  Moment 

Species 

p  X  10   c:r..^ 

€/k,  °K 

U,   Debye 

H^O 

2.75^^^ 

_ 

85.3("' 

1.81^ 

-r 

K 

2.6o(^) 

0.83p(^^ 

,,,W 

- 

oh" 

3.30^^) 

1.83(^) 

lilJ"'' 

- 

(a) 

Refe: 

rence  37        (b) 

Calculated  in 

this  work 

(c) 

Refe; 

rence  78 

TABLE  11 


PROPERTIES  AND  PARAMETERS  FOR  THE  SOLUTES 


Parameter 

Solute 

Hard  Sphere 

2.58 

He 

2.63 

»2 

2.87 

Ar 

3A0 

^2 

3.h6 

CH4 

3.82 

'h 

5-51 

neo-C^H.. 

l.kh 

Parameter 
e/k,    °K 

Pol. 
P  ^ 

ariz ability 

io2^..3 

0 

0 

6.03 

0.204 

29.2 

0.802 

122 

1.63 

118 

1-57 

137 

2.70 

200.9 

6.2l(^> 

232.5 


(b) 


10.36 


(b) 


(a)   Reference  80        (b)   Reference  8I 

All  other  parameters  have  been  obtained  from  reference  i^. 
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appropriate,  however,  to  take  o  as  twice  the  value  of  crystal  radius  of 
the  ion.   Since  crystal  radii  of  ions  are  difficult  to  determine  there 
is  considerable  disagreement  between  the  values  reported  by  various 
workers.   For  this  reason  a   values  for  K  and  OH  V7ere  obtained  as 
follows:   Equation  (173)  was  written  for  10  and  20^0  KOH  solution  at 
25  C  and  then,  by  using  a    as  calculated  earlier  and  0     as  2.58  A  , 

T 

these  two  equations  were  solved  simultaneously  to  obtain  o  tor  K  and 
OH  .   It  may  be  noted  that  these  equations  were  not  solved  in  a  rigorous 
manner  but  by  a  trial  and  error  procedure  using  o  for  ions  in  the  first 
trial  as  twice  the  value  of  crystal  radii.   The  values  obtained  are 
listed  in  Table  10.   As  in  the  case  of  water,  the  ionic  diameter  o 
obtained  by  the  above  procedure  varies  with  temperature.   These  varia- 
tions will  be  discussed  later.   The  values  of  ^  for  ions  as  obtained 
above  are  quite  different  from  hydration  radii  (Table  8) .   In  addition 
the  value  of  —  for  the  hydroxyl  ion,  I.30  A  is  much  larger  than  the 
value  G.hS   A  obtained  from  ionic  mobility  measurements.   The  abnormally 
low  value  of  the  radius  calculated  from  ionic  mobility  for  the  hydroxyl 
ion  is  due  to  the  special  mobility  ("proton- jump")  mechanism  of  these 
ions.   This  mechanism  was  discussed  by  Bernal  and  Fowler  (79)- 
c  for  Solute  Gases 

The  values  of  0  for  solute  gases  were  obtained  from  reference 
{kQ)    and  were  the  same  as  used  by  Pierotti  (37) .   They  had  been  cal- 
culated from  second  virial  coefficients,  and  are  tabulated  in  Table  11. 

Calculation  of  the  free  energy  of  introducing  the  solute  molecule 

c 

into  the  cavity,  g   ,  by  means  of  equation  (I52)  requires  a  knowledge 

CL 

of  the  energy  parameters,  e/k  for  all  the  species  involved. 
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c/k  for  Watv^r 

No  single  best  value  of  e/k  for  water  can  be  found  in  the 
literature.   Values  ranging  from  I67  K  to  77p  K  have  been  reported 
(37)  >  tiut  none  of  these  seem  consistent.   Pierotti  (37)  obtained  the 
value  of  e/k  for  v;ater  from  the  solubility  of  various  gases  in  water, 
using  equation  (128).   Since  we  are  using  the  same  general  theory  as 
Pierotti,  the  same  value  of  e/k  for  water  will  be  used.   This  value  is 
shown  in  Table  10. 
e/k  for  Ions 

The  parameter  e/k  for  ions  does  not  seem  to  have  been  reported 
in  the  literature.   However,  it  is  possible  to  estimate  this  parameter 
from  various  theories  of  dispersion  interactions  (^9,82).   The  theories 
of  London  and  Mavroyannis-Stephen  are  commonly  used  for  this  purpose. 
According  to  the  London  theory 


u 


,     .,  p  p.I  I .   , 

'^     _  -^  '  1  0,  1  /  1_ 


(  ^  )  (11^) 


aj    2  I  +  I.   •  o 
a    J    r 

where  u  .  is  the  dispersion  interaction  between  the  solute  molecule 
and  the  solvent  species  j.   I   and  I.  are  ionization  potentials  of 

solute  and  solvent  species,  respectively.   According  to  the  Mavroyannis- 
Stephen  theory 

1/2  2-  _ 

o.  J 

where  p   and  p .  are  respectively  the  polarizabilities  of  species  a  and 
j  in  the  mixture.   Z  is  the  total  number  of  electrons  in  a  particle, 
a_^  is  the  Bohr  radius  and  is  equal  to  O.5292  A  ,  and  e  is  electronic 
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charge.   The  dispersion  interaccion  from  the  Lennard-Jones  6-12  poten- 
tial is 

u^  =  -  U  .a°  .  (  ^  )  (176) 

Comparing  equations  (i?^)  and  (I70);  and  (175)  ^^i^  (176),  we  obtain 

he   .a°.  =  J  p  p.  _^-  ?  .,       (London)         (177) 
and  ^/2  2-  - 

6   ^  ^0    ^  ^0.^' 

^aj^aj   =  t  -f~U2 F~I7i  (Mavroyannis-Stephen) 

'^         j  (178) 

For  like-pair  interactions,  these  equations  become 

hco°  =   I  p^l  (London)  (179) 

keo°   =  7  aJ^Vp'^^^Z^^^       (Mavroyannis-Stephen)     (I80) 

To  compare  the  validity  of  the  above  dispersion  theories,  equations 
(179)  and  (16c)  were  used  to  calculate  e/k  parameter  for  molecules  such 
as  0^,  Ar,  ecc.   London's  theory  was  found  to  yield  values  which  were 
much  smaller  chan  those  commonly  reported  in  the  literature  (^).   On 
the  other  hand  the  values  calculated  by  the  Mavroyannis-Stephen  theory 
were  in  very  good  agreement  with  literature  values.   According  to 
Reed  (^9)j  fo^  monatomic  gases,  the  Mavroyannis-Stephen  expressions 
reproduce   the  empirical  parameter  e/k  quite  accurately.   Thus,  in  this 
v7ork  equation  (I80)  is  used  to  calculate  e  for  ions.   After  substitution 
of  the  values  of  a   and  e  this  equation  can  be  written  as: 
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3.11.6  xic-^V/^z^/2 

£  =  Z (181) 

a" 

To  calculate  e  for  K  and  OH  by  means  of  the  above  equation, 
V7e  need  the  polarizability ,  p  of  these  ions.   The  polarizabilities  of 
most  ions  have  been  tabulated  in  the  Landolt-Bornstein  Tables  (78). 
However,  the  polarizability  of  OH  was  not  available  in  the  literature, 
and  was  calculated  from  molar  refraction  data  by  the  following 
procedure  (83,  84): 

The  energy  quantity  corresponding  to  the  second-order  Stark 
effect  of  a  system  is  (83) 

AE  =  -  |pi'2, 
and  the  electric  moment  induced  in  the  system  is 

The  polarization,  P  in  unit  volume  can  be  expressed  in  terms  of  the 
index  of  refraction  n  as 

p.|.(^"-'4p£=fe4^£-  (182) 

'  n  +  2 

The  molar  refraction  R  is  defined  by  the  equation 

2       ,  ~ 

n  +  2    -* 
or 

p-St=2^='1°"''»-'  (183) 

Molar  refraction  data  for  OH  were  obtained  from  reference  (78).   The 
polarizability  value  of  OK  calculated  by  the  above  procedure  is  shown 
in  Table  10.   Also  included  in  this  table  are  e/k  parameters  for  K  and 
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OH  calculated  by  means  of  equation  (l8l). 
e/k  for  Solute  Gases 

For  most  of  the  solute  gases  e/k  parameters  were  obtained  from 
the  compilation  of  Hirschfelder ,  Curtiss  and  Bird  {k3) .      These  authors 
have  listed  Lennard-Jones  parameter  obtained  from  second  virial 
coefficient  and  viscosity  data.   The  former  were  used  in  this  vrark. 
These  parameters  together  with  the  polarizabilities  of  the  solute  gases 
have  been  tabulated  in  Table  11. 

_h        -^ 
To  examine  the  relative  magnitudes  of  g  ,,  „  and  g  /i^rp  of 

equation  (128)  and  to  compare  contributions  due  to  various  molecular 
interactions,  these  quantities  for  the  solubility  of  0   in  KOH  solu- 
tions are  tabulated  in  Table  12.   It  is  seen  the.:  the  free  energy  of 
making  the  cavity,  g  is  the  largest  positive  quantity  in  equation 
(128).   In  addition,  the  effect  on  g  of  adding  ions  to  the  system  is 

—S 

considerably  larger  than  the  effect  on  g  .   The  free  energy  of  intro- 
ducing  the  molecule  into  the  cavity  is  a  negative  quantity,  and  the 
largest  contribution  to  this  quantity  comes  from  nonpolar  interactions 
of  solute  and  water.   The  contribution  due  to  the  induced  dipole- 
dipole  interactions  are  relatively  small  for  various  potassium  hydroxide 
solutions.   It  may  be  noted  that  the  nonpolar  part  of  the  solute-ion 
interaction  is  quite  important  and  that  it  increases  rapidly  with 

potassium  hydroxide  concentrations.   The  variations  of  quantities 

-^,       -s 

g  /kT  and  g  /kT  with  potassium  hydroxide  concentration  are  shown  in 

Figure  2k.      The  absolute  value  of  both  of  these  quantities  is  seen  to 

increase  with  increase  in  potassium  hydroxide  concentration;  however. 
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KOH  Concentration,  Wt .  ^ 


J+0 


50 


Figure  2k. 


Various  Quantities  of  Equation  (128)  for 
Solubility  of  Oxygen  at  25°C. 
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—h  — s 

the  rate  of  increase  is  much  larger  for  g^ /kT  than  for  g  /kT.   Thus 

Cc  (X 

the  energy  Xi/hich  a  solute  molecule  requires  to  "dig"  a  cavity  in  the 
fluid  is  a  controlling  factor  in  determining  its  solubility.   Also  it 

is  seen  that  the  contribution  to  the  free  energy  which  allows  the 

RT 
solute  molecule  to  wander,  ln(  -rr  2  n.)  is  relatively  constant.   The 

.     \   y      j/ 

relative  magnitudes  of  various  quantities  of  equation  (128)  and  their 

concentration  dependence  for  other  gases  studied  in  this  work  were 

of  the  same  general  nature  as  for  oxygen.   It  may  be  pointed  out  that 

this  theory  is  capable  of  predicting  both  salting-out  and  salting-in 

effects.   The  latter  is  predicted  when  the  magnitude  of  the  free  energy 

of  introducing  the  molecule  into  the  cavity  is  more  than  the  sum  of 

the  free  energy  of  forming  the  cavity  and  the  mixing  term. 

The  values  of  InN  K  )  for  various  solutes  calculated  by  scaled 

particle  theory  are  compared  with  experimental  values  in  Table  13 . 

It  is  seen  that  at  25  C,  the  agreement  is  very  good  except  in  the  case 

of  very  large  molecules  like  SF^  and  neo-C_H   .   In  most  cases  the 

agreement  between  theoretical  and  the  experimental  values  of  the 

quantity  y   K  (inversely  proportional  to  solubility)  is  better  than 

15^.   In  the  case  of  S?^  and  neo-C  K   ,  it  should  be  noted  that  the 

o         p  J-^ 

theoretical  values  of  ln(7  K^ )  for  the  solubility  of  these  gases  in 
water  are  lower  and  higher,  respectively,  than  the  corresponding 
experimental  quantities.   This  trend  continues  even  for  potassium 
hydroxide  solutions.   These  discrepancies  can  be  explained  as  follows. 
Firstly,  since  the  Lennard-Jones  parameters,  a   and  e/k,  are  not  known 
accurately  for  various  substances,  they  might  be  significantly  in  error 
for  these  large  molecules.   Since  the  calculated  values  are  very 
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sensitive  to  Z   values,  a  small  error  in  this  parameter  can  result  in 
an  appreciable  error  in  ln(7  K  ) .   Secondly,  the  various  assumptions 
involved  in  the  theory  might  not  be  justified  for  these  systems. 
These  assumptions  are  discussed  below. 

It  should  be  recalled  that  the  only  assumptions  involved  in 

— h 
the  calculation  of  the  free  energy  of  making  the  cavity,  g  are  that 

(i)  the  molecules  are  composed  of  an  inner  hard  core  and  an  outer 

soft  part.   (ii)  The  hard  core  diameter  is  temperature  independent. 

By  calculating  the  solubility  of  a  hard  sphere  solute  of  diameter 

2.58  ■A  in  various  solvents  Pierotti  (36)  has  shown  that  these  assump- 

o 

tions  cause  little  error  at  25  C.   The  following  assumptions  were 

made  in  the  calculation  of  the  free  energy  of  introducing  the  molecule 

into  the  cavity,  g  . 

a 

(i)   The  pv_  contribution  to  the  free  energy  is  assumed 

_s 
negligible.   It  has  been  shown  by  Pierotti  (36)  that  v  is' a  small 

negative  number  for  weakly  interacting  solutes,  and  should  become 

larger  as  the  interaction  energy  increases .   Since  we  are  concerned 

— s 
vjith  a  pressure  of  the  order  of  one  atmosphere,  neglecting  pv  should 

cause  a  negligible  error. 

s 
(ii)  The  second  assumption  was  that  Ts   can  be  neglected  as 

s 

compared  to  ¥  .   It  is  difficult  to  determine  this  quantity,  but  it 
a 
s 

is  known  that  s  will  be  negative.   Since  the  system  becomes  more 

^  _s 

organized  if  the  interaction  energy  is  large,  Ts  will  make  a  signi- 
ficant contribution  for  strongly  interacting  solutes.   For  most  of  the 

solutes  considered,  the  interaction  energies  are  not  large  and  therefore 

s 
neglecting  Ts ^  should  not  cause  a  very  significant  error.   However,  for 
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SF^  and  neo-C  H,„  these  energies  are  quite  high  and  therefore  appreciable 

_s 

error  nay  arise  from  neglecting  Ts-  • 

(iii)   The  third  assumption  which  was  made  in  the  calculation 

of  g  was  that  the  radial  distribution  functions  g  .(r)  are  equal  to 
a  ccj 

unity  for  r  >  a   .   Since  g  .(r)  is  a  function  of  the  interaction 
energies,  this  assumption  should  not  be  very  good  for  strongly  interact- 
ing molecules  like  SF^  and  neo-C^H^^.   For  such  molecules  the  solvent 

a  ^    Id. 

species  cannot  be  uniformly  distributed  around  the  solute  and  some 

function  should  be  assumed  for  g  .(r). 

aj 

(iv)   The  fourth  factor  which  can  be  a  cause  of  discrepancies 
is  the  form  of  the  nonpolar  pair  potential  assumed.   The  Lennard-Jones 
potential  was  assumed  to  hold  for  all  nonpolar  interactions.   However, 
for  molecules  such  as  SF^  and  neopentane  the  pair  interaction  might  not 
be  adequately  described  by  this  potential. 
Predicted  Concentration  Dependence 

From  Table  I3 >  it  is  seen  that  the  concentration  dependence 
of  the  solubility  of  various  gases  in  potassium  hydroxide  solutions  is 
predicted  quite  accurately,  except  for  very  large  molecules  such  as 
neo-C  H.   The  main  factors  which  might  be  the  cause  of  the  discre- 
pancies in  the  case  of  large  molecules  have  already  been  discussed. 
In  all  cases,  both  the  predicted  and  experimental  values  of  ln(-y  K  )  in- 
crease  with  increase  in  potassium  hydroxide  concentration.   To  facili- 
tate the  discussion  of  the  concentration  dependeice  of  the  solubility, 
the  experimental  and  calculated  values  of  ln('Y^K  )  for  the  solubility 

'\Xf   tx 

of  oxygen  have  been  plotted  against  potassium  hydroxide  concentration 
(Figure  25) .   (The  agreement  between  theory  and  experiment  for  this 
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system  is  neither  the  best  nor  worst  for  the  system  studied.) 

Results  predicted  by  the  electrostatic  theories  of  Debye-McAulay 
and  Conway  are  also  included  in  this  figure.   It  is  seen  that  the 
scaled  particle  theory  gives  very  good  agreement  for  the  concentration 
dependence.   The  calculated  values  are  slightly  lower  than  the  experi- 
mental values  for  low  and  moderate  potassium  hydroxide  concentrations 
whereas  at  very  high  concentrations,  the  predicted  values  are  higher 
than  the  experimental  values.   The  discrepancies  in  the  results  can  be 
due  to  the  assumptions  discussed  earlier.   In  addition,  a  small  error 
in  the  choice  of  the  parameter  a   for  the  ions  will  introduce  a  large 
error  in  the  results  for  concentrated  solutions.   The  theory  of  Debye- 
McAulay  gives  much  lower  results  than  the  experimental  values.   On 
the  other  hand  the  theory  of  Conway  gives  fairly  good  results  up  to  20^ 
KOH  concentration.   At  higher  concentrations,  however,  the  expressions 

of  Conway  become  invalid  and  predict  a  negative  solubility 

S°  -  S 
(i.e.,  >  l) .   Discrepancies  in  the  results  calculated  by 

^2 
electrostatic  theories  were  discussed  in  section  5-2  in  terms  of  the 

assumptions  involved  in  these  theories.   The  concentration  dependence 

predicted  for  other  solutes  is  of  the  same  general  nature  as  for  oxygen. 
Predicted  Terr.perature  Dependence 

The  results  calculated  by  the  scaled  particle  theory  for  the 
solubility  of  various  gases  at  6C  C  are  shown  in  Table  13-   The  agree- 
ment between  the  experimental  and  the  calculated  values  is  quite  good 
for  most  of  the  solutes.   As  before,  the  agreement  is  poor  for  the 
case  of  SF^  and  neo-C_E   .   To  examine  the  temperature  dependence  of 

G  ^  Ltd. 
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the  solubility  results  for  oxygen  in  20^  KOH  have  been  plotted  against 
absolute  temperature  (Figure  26.)   It  is  seen  that  the  temperature 
dependence  predicted  is  too  large.   To  make  a  quantitative  assessment 
of  the  temperature  dependence,  partial  molal  heats  of  solutions,  /SB. 
were  calculated  by  using  relations  given  in  Chapter  2.   The  calculated 
and  the  experimental  values  of  AHp  for  0   are  tabulated  in  Table  Ik. 
It  is  seen  that  the  absolute  values  of  the  calculated  results  of  AH 
are  much  higher  than  the  experimental  values.   Similar  discrepancies 
are  present  in  case  of  other  solutes.   The  differences  in  calculated 
and  the  experimental  values  can  be  explained  as  follows: 

It  has  been  mentioned  earlier  that  the  hard  sphere  diameter 
for  various  species  is  a  function  of  temperature  (see  also  references 
(85)  and  (86)).   By  calculating  the  solubility  of  hard  spheres  of 
diameter  2.58  A  (zero  polarizability)  in  water  it  was  found  that  the 
theory  predicts  a  decrease  in  the  o  value  for  water  (see  reference  3?) • 

The  value  of  J  for  water  obtained  at  80  C  was  2.7^  A  as  against  2.75 

o 

at  25  C.   This  implies  that  a   for  all  the  species  involved  will  decrease 

slightly  at  80  C.   Although  this  decrease  is  very  small  it  results  a 
significant  error  in  the  results.   To  examine  the  effect  of  the  change 
in  a  values  at  higher  temperatures  for  various  species  the  solubility 
of  Op  in  20^  KOH  solution  was  calculated  by  allowing  o  to  vary  with 
temperature  (see  Table  I5).   The  80  C  values  shown  for  a  for  ions  and 
oxygen  vjere  estimated  by  assuming  their  diameters  to  decrease  with 
temperature  in  a  proportionately  similar  way  to  that  of  the  water 
molecule  diameter.   Appropriate  values  of  G  for  temperatures  intermediate 
between  25  and  80  C  were  obtained  by  interpolating  between  the  values 
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Figure  26.   Effect  of  Temperature  on  Activity  Coeffi- 
cient of  Oxygen  in  20/o  KOH  Solution. 
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TABLE  Ik 
/^H  FOR  0/20^  KOH  SOLUTION,  cal,(g.mole)"''" 


Temperature  Expt .  Theory 

o„ 


25 
40 
60 
80 


-1100 

-2020 

-  812 

-1770 

-  655 

-1325 

-  kk6 

-1097 

TABLE    15 
PAR.^METER  a  FOR  SOLVEOT   SPECIES,    A° 


Solvent  Species 

25°C 

°2 

3.k6 

H^O 

2.75 

2.60 

oh" 

3.30 

80°C 


3.i^5 
2.7^ 

2.59 
3.28 
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shown. 

The  results  obtained  are  plotted  in  Figure  26.   It  is  seen  that 
these  a  values  lead  to  a  more  correct  terr.perature  dependence  for  the 
oxygen  solubility.   In  the  case  of  SF^  and  neo-C_H, ^,  the  decrease  in 

o  5  •'-'^ 

G  will  be  still  larger  and  will  therefore  introduce  appreciable  errors 
in  the  results  at  higher  temperatures,  as  is  observed  in  Table  13 . 

3 -U Comparison  of  the  Electrostatic  Theories  and  the  Proposed 

Cavity  Modal  Based  on  Scaled  Particle  Theory 

(1)  The  most  striking  difference  between  the  two  principal 
types  of  theories  is  that  whereas  the  Debye-McAulay  and  Conway  theories 
can  satisfactorily  predict  salting-out  effects  only  up  to  5^  ^nd  20^ 
KOH,  respectively,  the  proposed  model  gives  very  satisfactory  results 
for  potassium  hydroxide  concentrations  up  to  5C/o. 

(2)  The  electrostatic  theories  enable  one  to  calculate  the 
solubility  of  the  solute  in  an  electrolytic  solution  only  if  the 
corresponding  solubility  in  water  is  known.   On  the  other  hand,  the 
proposed  cavity  model  gives  the  former  quantity  directly. 

(3)  The  tenperature  dependence  of  the  solubility  can  be  fairly 
satisfactorily  calculated  by  the  cavity  model  whereas  electrostatic 
theories  cannot  attempt  such  a  calculation  because  the  parameters 
needed  for  these  theories  are  usually  not  available  at  higher  tempera- 
tures . 

{h)        The  proposed  cavity  model  should  be  able  to  predict  the 
observed  salting-in  phenomena  in  cases  when  the  total  molecular 
polarization  of  the  solute  is  less  than  that  of  the  solvent.   In  such 
cases  the  electrostatic  theories  predict  that  salting-out  should  occur. 
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In  fact  Bockris  et_  al.  (l6)  recognized  the  importance  of  nonpolar 
interac;;ion£  in  a  so -_:e /electrolyte  system  and  accounted  for  the  salt- 
ing-in  of  benzoic  acid  by  tetra  alkyl  animonium  iodides  in  aqueous  and 
ethylene  glycol  solutions  by  considering  Van  der  Waal's  forces  between 
the  solute  and  the  ions  in  addition  to  the  coulorabic  forces . 

(5)   The  electrostatic  theories  are  difficult  to  use  because 
the  parameters  involved  are  difficult  to  obtain  accurately.   On  the 
other  hand,  the  cavity  model  contains  parameter  which  are  either 
easily  av-iiable  or  can  be  calculated  with  sufficient  accuracy  and 
withouc  undue  difficulty. 

In  addition  to  the  above  obvious  advantages  of  the  cavity 
model,  based  on  the  scaled  particle  theory,  over  electrostatic  theories, 
there  is  an  important  conceptual  advantage:   The  proposed  model  is 
simple,  easy  to  visualize  and  starts  with  fairly  rigorous  equations  of 
statistical  mechanics.   The  solubility  of  a  gaseous  solute  in  organic 
solvents  (36),  water  (5?)  and  electrolytic  solutions  can  be  explained 
in  terms  of  single  theory  involving  no  assumptions  concerning  the 
structure  of  the  solvent.   Electrostatic  theories,  on  the  other  hand, 
start  wich  the  basically  incorrect  assumption  of  a  continuous  medium, 
are  difficult  to  conceive  and  becon.a  still  more  obscure  when  the  ideas 
concerning  the  structure  of  the  solvent,  hydration  shells  and  dielectric 
saturau^on  effect  are  introauced. 


CHAPTER  6 
CONCLUSIONS  AND  RECOMMENDATIONS 

An  experimental  study  of  the  solubility  of  nonpolar  gases  in 
potassium  hydroxide  solutions  has  been  made.   The  variables  considered 
are  the  various  physical  characteristics  of  the  solute  gas,  electrolyte 
concentration  and  temperature.   This  study  has  been  used  to  test  the 
electrostatic  theories  of  salting  phenomena  and  the  proposed  cavity 
model  based  on  scaled  particle  theory.   The  proposed  model  offers 
several  advantages  over  existing  theories  of  gas  solubility  in  electro- 
lytic solutions  and  is  shown  to  predict  fairly  accurate  values  of  the 
solubility  of  nonpolar  gases  up  to  an  electrolyte  concentration  of 
50/^  by  v;eight  and  for  temperatures  up  to  80  C.   The  proposed  model 
raises  an  interesting  point.   The  fact  that  the  ions  are  charged  does 
not  appear  to  affect  the  solubility  process  to  a  significant  extent. 
The  difference  between  the  solubility  phenomena  in  electrolytic  solu- 
tions and  pure  water  lies  in  the  fact  that  the  former  contains  a 
number  of  particles  (ions)  having  molecular  parameters  which  are  differ- 
ent than  those  of  water.   The  nonpolar  interactions  between  the  solute 
and  the  ions  are  found  to  be  the  only  interactions  which  are  important. 
This  is  in  contradiction  to  earlier  viewpoints  (I9)  in  which  the 
electrostatic  interactions  and  hydration  effects  were  supposed  to  be 
mainly  responsible  for  salting  phenomena. 

For  an  exhaustive  test  of  the  proposed  cavity  model,  systematic 
experimental  studies  of  the  kind  done  in  this  work,  should  be  made 
for  both  polar  and  nonpolar  gases  in  several  other  electrolytes.   The 
electrolytes  should  be  chosen  in  such  a  way  that  the  ionic  species 
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involved  have  widely  varying  molecular  parameters  and  valence.   Similar 
studies  should  also  be  made  to  observe  salting-in  phenomena.   It  will 
be  interesting  to  see  if  the  theory  can  accurately  predict  solubilities 
as  a  function  of  the  type  and  valence  of  the  ion.   A  severe  test  of 
the  theory  will  be  the  prediction  of  salting-in  effects.   The  theory 
might  also  be  used  to  predict  the  solubility  of  gases  in  molten  salts 
and  liquid  metals.   Finally,  the  thermodynamic  properties  of  ions  in 
water  or  other  solvents  should  be  amenable  to  study  by  including  terras 
involving  ion-ion  and  ion-dipole  interactions. 

The  theory  is  by  no  means  complete  and  leaves  scope  for  modi- 
fication.  Suitable  methods  of  obtaining  accurate  ionic  diameters  should 
be  developed  and  should  be  used  in  this  theory.   In  this  work  these 
parameters  have  been  calculated  from  the  solubility  data  themselves.  For 
solute  gases,  it  would  be  preferable  to  use  parameters  J  and  e  which 
had  been  smoothed  before  use  (o7) •   This  may  improve  the  predicted 
results  for  large  solutes.   To  improve  the  theoretical  temperature 
dependence  of  solutility,  a  reliable  method  of  predicting  the  tempera- 
ture dependence  of  C  should  be  developed.   A  further  refinement  would 
be  to  find  a  simple  expression  for  the  radial  distribution  function 
for  gas/electrolyte  systems,  and  this  function  should  be  used  to 
calculate  the  solubility  results.   Finally  a  means  should  be  devised 
to  take  into  account  the  entropy  associated  with  the  charging  of  the 
soft  potential. 
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APPENDIX  1 
SURFACE  EFFECT  AND  SOLUBILITY 
Intuitively  it  is  seen  that  the  pressure  inside  a  gas  bubble 
is  higher  than  that  over  a  plane  surface;  i.e.,  the  chemical  potential 
inside  such  a  bubble  is  higher  than  that  over  a  plane  surface.   Quanti- 
tatively, the  equation  for  the  energy  change  is  (88) 

dG  =  -SdT  +  VdP  +  adA  +  Z  jLtr^  dn.  (iQk) 

The  above  equation  is  valid  for  curved  surfaces  under  the  assumptions 
(i)  the  surface  tension  is  unchanged  by  the  curvature   (ii)  the  inter- 
face within  the  system  does  not  change  the  pressure  on  the  exterior. 


ial  ^(^) 
pertains  to  transfer  of  matter  when  the  surface  area  is  unchanged,  i.e.. 


It  may  be  noted  that  the  chemical  potential  ^.    in  the  above  equation 


iS.^^  =    (  S~  I  ^  ,  „  (185) 


'  i  i   T,P,A  ,n. 

s   J 

For  a  spherical  bubble,  however,  the  area  will  change  according  to 

the  transfer  of  matter.   Thus  the  volume  change  or  addition  of  dn.  moles 

of  various  components  is: 


dV  =  Z  V . dn . 

X   1 


For  a  spherical  bubble: 

dV 

dA  =  2  -^ 

s      r 


V. 
=  2  Z—  dn.  (186) 


r   2. 
Thus,  for  transfer  from  a  curved  surface,  we  can  write 


ito 


I4l 


dG  =  -SdT  +  VdP  -^  c(2  2—  dn.)  +  2  M.  '^  dn. 


(p) 


r   1 


=  -SdT  +  VdP  +  Z  {jj}.^   +  2  -^  )dn_.  (iS?) 


1 
1 

We  note,  therefore,  tha^t  the  chemical  potential  for  a  curved  surface  is 

V.a 
raised  by  an  amount  2  over  the  plane  surface  chemical  potential, 

(P) 

u.    due  to  surrace  errects. 

(c) 

The  chemical  potential  for  a  curved  surface  jj..        is  given  by: 

Mi^>=.r'*2^'  (186) 


By  equation  (S)  of  Chapter  2 


f .      -^    o 
In 


(P)    RT   RT 

and 


7.x. 


(p)     ^i        r  u°  -  u""  -]  ,.,.. 

Now,  in  the  case  of  surface  effects,  equation  (l)  of  Chapter  2  and 
equation  ( l58)  give 

G    (C)    (P)   ^  \^ 

V.a 
=  ^°  +  RTlnf.  -r   2  -^  (IQO) 

The  equation  similar  to  (I89)  for  a  curved  surface  can  be  written  as 

/p\   f .     -       .        V.a  _ 
^^Cj  ^  ^  ^^p  I  (^o  _  ^  ^/^^  ^2^1  (191) 

(C) 

where  x.    is  mole  fraction  in  the  case  of  a  curved  surface.   From 
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equations  ( I89)  and  (I9I)  we  have,  then 


1  II 


e^P  !  7^^  I  (192) 


(P)  -  "="^  ;  rRT 
X.         *- 


Sample  Calculation: 


0^/Water  System 


Cf  =  79 -7  dynes /cm.  for  water 


Vq  =  32.00  cn^yg.mole. 

2 


=  exp  (^  — ^ j 


X 


r  =  0.1  cm . 
(C) 

K^  =  1.C007 

X 

%   Error  =  0.07^ 

Radius  of  bubble  ^  Error  in  solubility 

r  ,cm  . 

0.1  0.07 

0.01  0.49 

0.001  k.66 


Rep 
10 


APPENDIX  2 

ESTIMATE  OF  ACCURACY  OF  EXPERIMENTAL  SOLUBILITY 

MEASURE>2XTS 


proQucibility  in  the  concentration  range  of 

'-'   to   10"'   g.mole  of  solute   gas/litre  =        ±1.05o 


Reproducibility  in  the  concentration  range  of 

10"5  to  10"°  g.mole  of  solute  gas/litre  =  ±5.0fo 

Chromatographic  error  =  ±0.2^ 

Saturation  error  =  ±0.2^ 

Injection  error  =  ±0.25& 

Error  due  to  KOH  concentration  measurements  =  ±0.1^ 

Miscellaneous  =  ±0.3^ 

Total  accuracy  for  dilute  solutions  =  ±2.0^ 

Total  accuracy  for  concentrated  solutions  =  ±6.0^ 


This  is  the  standard  deviation  from  the  arithmetic  mean 


li^3 


APPENDIX  3 


^   ,   l^    ,   l\S     FOR  HYDROGEa,  HELIIM,  ARGON, 
SULFUR  HEXAFLUORIDE,  METH.ANE,  AND  NEOPENTANE 
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Figure  27-   Partial  Xolal  Heats  of  Solution  for  Hydrogen. 
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Figure  29.   Partial  Molal  Heats  of  Solution  for  Argon. 
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Temperature ,   K 

Partial  Molal  Heats  of  Solution  for  Sulfur 
Hexaf luoride. 
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Partial  Molal  Heats  of  Solution  for 
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Figure  33-   Partial  Molal  Energies  of  Solution  for 
Hydrogen. 
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Figure  36-   Partial  Molal  Energies  of  Solution  for 
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